CHAPTER XXIV. 


EULERIAN INTEGRALS, GAUSS' II FUNCTION, ETC. 


853. The Original Forms of the Eulerian Integrals. 
The properties of the two important integrals 


„рє; 
-lda 1 1\2 
= (8) | апа = |2) | (log 7) dz 
(1—2%) * 
were the subject of several remarkable memoirs by Euler. 
His investigations were published in the Jnstitutiones Calculi 
Integralis, 1768-1770, and are of great importance in the 
general theory of Definite Integrals. The notation above, viz. 


(P and Hi is that used by Euler, and the above forms are 


those in which the integrals were studied both by Euler 
and Lagrange. In each of these the value of the integral 
was supposed to change by the variation of p and q; the 
^ which occurs in the first integral was supposed to be a 
constant. 

Legendre, for the purpose of charaeterising these integrals 
and honouring their great discoverer, named them “ Intégrales 
Eulériennes."* The second part of Legendre's Exercices de 
Calcul Intégral is devoted to a discussion of their properties. 
He adheres to the notation (8) for the first integral, but 
suggests the notation г(?) for the second, regarding Г(а) as 
a continuous function of a. 


* Exercices de Calcul Intégral, par A. M. Legendre, 1811, р. 211. 
49 
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854. The More Convenient Modern Forms. 

The above forms of the integrals are not the most convenient 
in practice. Taking the first integral, write z^—y, and put 
p=nl, q—nm. 


Then 
2-1 12 
n Жу, 
1 2215 19 Yin 
0—2) * ба“ ^ 


Jo 
Taking the second integral and writing log += y, that is 


x=, and putting : =n, 


Ва 


Hei (log 2 ds | езуд. 


0 

855. Definition. 

We shall therefore define the FIRST AND SECOND EULERIAN 
INTEGRALS as 1 
B(l, ")= | g!-1(1—2)n-1dz 

0 
and г) = | е-*т"-1йх, 
0 
and refer to them respectively as the BETA and GAMMA 
Functions. This is now the commonly accepted notation and 
nomenclature. 

856. In Gregory's-.Examples (p. 470), the digamma F(l, т) 
is used to denote what we have above defined as the Beta 
funetion. It will be observed that B(l т) is m times the 


integral discussed by Euler, that is (P) 


We shall assume in our subsequent work that all the quanti- 
ties l, m, n are positive but not necessarily integral, and 
further that they are real unless the contrary be expressly 
stated. 


857. The Beta Function is symmetric in / and m, that is, 
В(/, m) = Bm, I). 
If in the Beta function 
1 
B(, т)= | 2-1 (1 —2)"-14х 
0 
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THE BETA FUNCTION. 51 
we write 1— y for z, we obtain 
Bi, m=- (1—уу-!у"-1йу= [| m0 ууу 
=f g"-1(1—2y-!dz-—B(m,l) ; 


whence it appears that B(l т) is a symmetric function of 
l and m, the l and т being interchangeable and 


В т) = B (m, 1). 
This property might be exhibited by writing B(l, m) as 


ВЦ, т)=» | [9-1 (1 —2)n-1-- 2-10 —2) 1]dz. 


858. Case when / or m is a Positive Integer. 

When either of the two quantities $, m is a positive integer, 
the integration is expressible in finite terms. 

Suppose m is a positive integer, 


1 
BU, m)=| gl 1 (1 — a yn-1dz, 
0 
and by continued integration by parts 
2! gi 
-|7 (1 ay ey M- а в)" 


qii? 
+ TTI U 1)(т—2)(1—2)”-3-+... 

gitm i 
THE ra" De 2.2.1] 
(m— 1)! 


~ 1041)... (т 1): 


Similarly, if | be a positive integer, 


Bi. m) (1—1)! 
miser) ... (m+1—1)’ 
and if both be positive integers, 
_ (1—1)! (m— 1)! 
ть Г: 
859. Various Forms of the Beta Function. 
The Beta function may be thrown into many other forms 
by a change of the variable, and therefore many other integrals 
are expressible in terms of the Beta function. 
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Thus: (1) Let y=. 
Then B(l, m) f y= (1—y)n-1dy 
0 
a т 1—1 £ m—1] 
=|, Q-» ie 
=. 1 zi gl-1 m-1q. 
=з]. (а—х) g. 
Hence f gl-!(a—z)"-1dz—a*"—-'B(l, m). 
0 
2) Let Me. 
(2) e. e 
1 
Then Bil, m)=| y (1 — y) My 
0 


© 


1 йы da: 
-f аа) hara 
gm-i 

T o (l+a)+™ А 


and since l, m are interchangeable this must also 


0 11-1 
sharr” 
which would have appeared immediately if we had made the 
e EU 1 
substitution к= instead of fe 
Note also that the symmetry in l, m may be exhibited as 
| w LE gi-1--3m-1 i 
Bil, m)25] ИТАК”? 
whilst for all positive values of Г and m we have 
1 aj—1. gmt 
|, ттен» 0 
So that, for instance, 
z(1-2*).,. o. 128 (0 4 x), 
|, caps 4—0: and Í ( rege d 2B (6, 12), 
(3) Putting 


x 


y up pi. ах 
1+9 a+a’ WET) (= а)’ 
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1 
Bü, m)- | уу)" 
1 т 1—1 | ee т—1 ах 
нея 1—1 т— BN oot н, 
= |, (1+9) E p care, ris (ха) 


1 gi-1 ]— m-—1 
атау || рн dr 


1g-(1—29-*. BG, 9) 
Hence | ео атаў 
This is Abel’s transformation ((Zuvres; Vol. I., p. 93). 
2—0 
(4) Put y=. 


Then B(l, т) =| у—1(1—у)"—14у 
А ON гаа 
= рр |, (2—5): -1(a—2)"-!dz, 

апа | (ж—Ь)!—1(@—хж)"—14ж=(а—Ь)+"-—\ B(L т). 

Неге the limits have been changed to any arbitrary con- 

stants a and b. 

(5) Transform by the formula 954—0 
Here the limits remain unaltered, for if y=1 we have z— 1, 
and if y=0, х=0. 
ВО, m) [ya — dy 


a f { ba di { | a(1—2a) " ab dx 
o Ve 4- (b— a )a: a+(b—a)a {a+(b—a)x} 
и 1 ge (1— zy -1 
p. v). Фа o^ 
1 йт 1(1—ж)"- -1 
Ben Grea cm 
E т -1(1—ж)"- 1 
also obviously | op eva Еа Бау = 
and if we write a—b=c, 


141-1 (1 =g)” 
у толя 


В( т); 


Bi, m); 


В(/, т). 
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54 CHAPTER XXIV. 
(6) In the last transformation, put z—sin? Ө. 


чы 2 6119-20 cos?™-2 0 i 
Then ji (a cost OLD aint OF 2 sin Ө cos 0 49= p BC, m), 


; [ sin*-! 0 сов*"-1 0 de— 
p (а cos? 0-4- b sin? Өу+" 2 aD 


l, m, a and b being positive constants. 


В( m), 


т 


z 
(7) I =| віп? 0 соз? Ө 00 is expressible in the same way in 
0 


terms of a Beta function. 


Let sin@=J/z, ie. cos d0— Lae Ч 
9 Ja 
rafa a? (1— ar n 
0 Wr 


рт. ү 9+1 | 


=5[ 2° (1—2) * ах 
0 

_1р(р+1 9+1 
-gb( go y) 


This also follows from No. (6) by putting a=b=1. 


860. Properties of the Gamma Function. 
Consider next the Gamma function, viz. 


T(n) =f EERE: 
* Integrating by parts 
Е (n)=| ие +(n— nf enondo; 
- 0 0 


and whatever n may be, provided it be finite and >1, 
—g"-1e-* vanishes at both limits. 


Hence Tl'(n)—-(n—1)I(»—1) 
Similarly, (п —1)2(n—2) T (n — 2), 
and so on. 


in the ease then, where n is a positive integer, 


I'(n)2(n—1)(n—2)(n—3)...3.2.1T'(1), 
and га) e &=| -e| =1, 


whence Г (2) =(2 — 1)! in that case. 
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861. Working Properties. 
We then have the properties 


Ты вороны: нь... Т. 
Dd dat i. заочно dois i II. 
and when т is a positive integer, 
Dh dmm Jeune nnm ПІ. 
The Gamma functions of the positive integers are then 
Г(1)=1, 
Г(2)=1.1=1, 


Г(3)=2Г(2)=1.2, 

Г(4+)=3Г(3)=1.2.3, 

Г(5)=4Г(4)=1.2.3..4, 

etc., 

from which a rough idea of the march of T (Œ) as a continuous 
function may be inferred, viz. à minimum existing somewhere 
between х=1 and 2=2, and then after z—2 a quantity 
increasing more and more rapidly. 


862. In any case the equation I'(n4- 1) =т Г(т) furnishes a 
means of reduction of the Gamma function of any number 
greater than unity to а Gamma function of а number less 
than unity. 


For instance 
Гоу) e Гмм) =. V Го). VE TG) o A398: 8T) 
z1í,.11,2.5.8 T'(3). 

That is, the Gamma function of any number greater than 
unity can be connected with the Gamma function of а number 
which is not greater than unity ; so that it is already obvious 
that when we come to the caleulation and tabulation of the 
numerical values of Gamma functions it will be unnecessary 
to tabulate T(x) for any values of ж except those which lie 
between 0 and 1. 


863. A Caution. 
The student should guard against the idea that the equations 


re e-*vw-!dv and I(¢#+1)=aI (ze) 


are co-equivalent. They are not so. The latter is а conse- 
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quence of the former, not the former of the latter. The latter 
is а functional or difference equation, viz. 

$(xr-1)—z$9(x) or wuw,4-2zu;, 
and such equations may have many solutions. What is proved 
is that w=] e 4*-! dv is a particular solution оў u y; = 201. 


0 © 
But so also are 4| e-*V*-dv when А is any constant, or 
0 


such an expression as 
A+B cos* 27a f” 
C+D sin$ m 
where А, В, C, D are constants, for these multipliers are not 
altered when х is increased by unity. Мог does it follow 


e~ 2—1 dy 


that | е—° 1—1 dv occurs as a factor in all solutions of the 
0 


difference equation. 
The solution of 4,,,—2w, is obviously 
Áz(z—1)(x—2)... (r3- 1)ru, 

when А is either a constant or some arbitrary periodic function 
of x whose periodicity is unity, and which therefore does not 
alter when 2 is increased or decreased by any integer, and u, 
any assumed initial value of и. We shall return to this 
matter later. 


864. Transformation of the Gamma Function. 
As in the case of {Һе Beta function, transformations of the 
variable will give rise to other integrals. 


(1) We have seen that z—log: ог у=е-* produces 


oo 1 1\"-1 
Г(п =| ж"-—1 da log — dy, 
бу] arte (бов) dy 
the form studied by Euler. 
(2) If we write kx for 2, 


г(ш)= | ееп dg ; 

0 

whence [em gh-l d; L9 а) 
x 0 v 


provided k be а real constant (see Arts. 1159 to 1162 
and 1327). 
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(3) If we put z^—7 where n is positive, 
год, enar: 
es | ow ду=т Г(п)=Г(®- 1). 
0 
In this case, if we put п=$, 


0 0 
and this leads to an easy calculation of T (4). 


For ü г(р#=] е-”йвх f е-”ду, 


and аз z and y are independent variables and the limits 
constant, we may write this as 


Fee ад). 
0*0 


Now, regarding z, y as the Cartesian coordinates of a point 
we have to sum all such elements as e- @ +”) да dy through an 
infinite square in the positive quadrant, two of whose sides are 
the coordinate axes. 


J 


о x 
Fig. 313. 


Transforming to polars, we have to sum 


ест 60 дт 
through the same square. 
Let z—a, y=a, where а= о , be the other two sides of the 
square. Then for the portion of the square which lies inside 


the circle z*-- у2= а? the limits for Ө are 0 and 5, and for т 
0 and œ. 
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Hence the portion within the circular quadrant contributes 


TET s di UD. E T ОРЕ АЙ, 
NE rdrdó- | ve dr=7| ie | ne a 


At points of the square outside the circle the elements are 
never greater than e-^'r 60 ôr, and when о is made sufficiently 
great this becomes an infinitesimal of higher degree than the 
second, and hence in the double integration disappears. There- 
fore the portion of the area between the circle and the square, 
exterior to the circle, contributes nothing. 

Hence the value of Г($) is + Ут, and as all the Gamma 
functions are from the definition essentially positive quantities, 

Г()=Ут.* е 
865. We may also regard the investigation of | e~“ du as 


the problem of finding the volumet bounded by the plane of 
x-y and the surface formed by the revolution about the z-axis 
of the curve 2=е-*', for this volume may be regarded as 


Fig. 314. 


being built up of cylindrical shells whose axes coincide with 
thez-axis. "The volume of this solid is then | 2rudu.z, where 
u is the radius of a section parallel to the т-у plane, 


0 
2» | тет du. 
0 


* Euler, Tom. V., des anciens Mémoires de Pétersbourg, p. 44. 
T Airy, Errors of Observation, р. 12. 
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But dividing it by planes parallel to the coordinate planes 
of z—0 and y=0, the volume is also expressed by 


FA eem] а= oe dex | еду 
-([ 7 dz); 


whence | е- dr E 
0 2 


This gives another geometrical interpretation to the work 
of the preceding article. 


866. When n is diminished without limit | e-*a"-1 dæ be- 
0 
comes infinite. For the formula T (n+1)=nT (n) holds for 
all positive values of n. Hence 
Lt, I (n) м 0, 
ie. T(0)—oc. 
This is also obvious from the integral itself. For the 
integrand - (for the сазе n=0) takes an oo value at the 


lower limit, and the principal value of the integral becomes 
infinite (see Art. 348). 


867. Connection of the Two Functions. 
We shall next prove that the Beta function is expressible in 
‚ terms of Gamma functions, the connection being 


ГО). (m) 
BC, ПЕ ; 
Consider the double integral 
IUE e7 (ay) -! хе-=хтах dy 
o Јо 


[that is wy is written for 2 in the integrand of I'(D), and this is 
multiplied by the factors of the integrand of T(m+1)], хе 


r=" Г сит) +т- цу -14у ах. 


Integrating first with regard to 2, we have 


Til+m) 

I-| y d 

|, o Eads 

=I (l+m)B(l, m), by Art. 859 (2). 
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But changing the order of integration, taking y first, 
I =f | e-ayltm- TORO dy 
=f ganent i Oa, 
0 а? 
=T of елш 
=T (l) К (m). 


Hence B(} му L0 Lin 


T4m)' 
868. Deductions. 
It further follows that 


_ I'(Id- m) T(n) 
Bin 2» туй 
and therefore that 


Bern) рна) dy EAE CREE ES) 


I'l4-m-r0) ' 
which is a symmetric function of /, т, п. Hence we have 
В( т) B(l-+m, п)=В(т, п) B(m+n, 1) «B (n, l) B(n-F-L т). 
Hence also 
BL, m) Вт, т) B(L4-m 4-0, py FUTT C XD), ete 
869. It now follows that the results of the transformations 
of the Beta function given in Art. 859 could be further expressed 
ағ Gamma functions. 
Thus 
(eee = 1 Bü m= 1 T'() P (m) 
о (b-+-ca)it™ stica E btob T (L4-m) 


f sin?/-19 cos?n-19 49 Ретин іа 1 Г(0Г(т) 
o (a cos? 0 +b віп0) + Ep а"! T (L4-m) ' 
күзүк 
удай ODIT 
Li 
| sin? 0 cos? 40-5 В(Р+', in HG) 
| ТР вены) 


ete. 


The last of these integrals has already been used in earlier 
chapters, for convenience of calculation, with a temporary and 
limited definition of Г. 
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"жава 4 


THE INTEGRAL ГЕ "s 
870. We have also in Art. 859, Case 2, the integral 


A. pe T(z) P (m) 
‚ rage d BU т) Tm)’ 
Then, since Г (1)=1, we have 

-1 


Put [+ т=1. 
11.9, 


Рога Lu a 


where m is a positive proper fraction 
We have then to consider this integral next. 


871. The Integral 7 — i € where 0 «n <1. 
The integration | may be separated into two parts, viz 


| 


In the second part put BE 


Then 
oq" са xil. йы hh aan 
аа de itr )ay= =| E ga d 
y 
Ma i. a 
Hence I zi ifs da, 
and by division 
i=l- — 1422—13... ( —1)*--(— perf 
Hence 
I= (а-а) ея Oc (— 1)5) 
Tu х" -14 g7” 
aa aoe) 
] 1 1 1 
lr LER state м n k+n 
1 1 1 1 
Наив ee ay 
1-х? 
— 1)% а ed 
ceo intus 
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Now cosec z= 
1 1 1 1 1 1 1 
БОКЕ 3:4 Ded рва UE 


(Hobson, Trigonometry, p. 335.) 


Hence 
1 1 1 1 1 1 1 
a; CTR т-а Ше Саца са 0" 
— т е 
віп тт’ 


and since in the limit when / is made PR Y large the 


last term of the series for J, viz. (— 1)* БЕТТ becomes zero, 


at 


the portion of I within the brackets becomes —— — 
sin nr 


grit an 

1+2 
note that аз z lies between 0 and 1 and is а positive proper 
fraction, z**! is diminished indefinitely by an infinite increase 
in k. If then this integration be expressed as a summation 
according to the definition of Art. 11, each term of the sum- 
mation is diminished without limit, and may be regarded as 
an infinitesimal of the second or higher order when Ё is 
sufficiently increased. 


1 g^-1J-g7^ 
И ae ^m 
Непсе Linef 2 а dz=0, 


Also as to the remainder, viz. f gett dz, we may 
0 


and we are left with 
SÉ ras (№ 
от“ sinn 
872. An Important Result. 
It now foliows that 


where 0<n<l. 


l(»)l'(1—2)— “еее (0<n<1). 
As a particular case put n=}. 
ГАТ Ay я 
(r(5)} uem ији, 
sin 5 
2 
and T (3) — Jr, аз has been seen before, Art. 864. 
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Again, put n=}. 


r(Dr(--—-«3. .. ЕЕ 


: "a 1 г(2) 
Put n=}. 
г(в)г( = =з RUN т etc. 
o- 


Hence Г ($), Г($), etc., are expressed in terms of Gamma 
functions of numbers which are <4; whence it will appear 
that if all Gamma functions were tabulated from T (0) to T (4), 
all others could be found by this theorem, together with the 
theorem Г (n+1)=ĦnT (n). 

The result Г(®)Г(1—”)=- еј 

sin т 
in an earlier chapter, Art. 592, in the calculation of a certain 
arc of а Lemniscate. 


Since ы Кы nl (n) and Г(п)Г(1—»)=_——, 


Sin 17 


‚ was temporarily borrowed 


this formula may be written 


T(124-2) Г(1— т) =; 


Sin NT 


Arr Were, шефу 


873. To show that 


Ше =” сы 


We are now able to consider the continued product 


1 2\ ./3 n—1 
Per GGG Pa) 
where n for the present is any positive integer. 
By writing it down again in tke reverse order, multiplying 


the results, and noting that 


O CA P («n 


n n 


F 
Sin — т 
n 
т т т т 
we have : P= ^ —— 9 — ——— à 
ке MEME SEUEUINTUM uiu 
sin я sin * В | sn —— 
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and since 
sinnÓ и: т\з 2т\ . 37 Е 
roa lsin (0+2) sin (9+ ==) sin (e) ... Sin (0+ 
(Hobson, Trigonometry, р. 117), 
we have in the limit when 0=0, 


n—1 T) 
n 


l8 93:5. Эт . (п—1)т 
n=2"-1 sin — sin — sin — ... sin (ir 
n n n n 


т 


n—1 
2 27-1, and P being positive, we have 


Hence P?— 


n-i 


т) 9) т). 


ni 
874. Gauss’ IT Function. 
Taking the original Eulerian form of the Gamma function, viz. 


T(n) af (log D) de, 


x 


E 
1—2* 


and remembering that Гі, -„ =log * (Diff. Calc., Art. 21) 


Е 
È и 
we may write 
iex fae) 
(ов) = TUA ы 
„ЩЩ, 
where e is somethiüg which vanishes in the limit when и 
becomes infinite. 
Let us take и as а positive integer. 


1 1 п 1 
Then го) | ui (1-9) ic | є dx. 
0 0 
In the first integral put z—y". 
1 1 
Then l(n)— “| y0—-yyàye| e ах, 
0 0 
and as y is a positive integer, 


E —1 n-lg,— ( —1)! 
| y" (1-9) LE E i P ҮЗА, (Art. 858). 


n (uI)! | 
Hence T'(n)=p nar 
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Hence, making u increase without limit, the integral ulti- 
mately vanishes, and 
ely 
Pes Door tek 
(n)= Lt, e i п(п- 1)... (n3-4—1) 
or, which is the same thing, 


" A u! i 
T(n) = Lt, - 04 n(n4-1) ... (n3- y. —1)' 
and writing n+1 for n, 
П Рр oo 44 


Г 1)= Lt, -o u” : 
(ma od e" ede Dor (ib) 
This limit is known as Gauss' II function, and is written 
1:2:8 ip 


TE Meme ET) (n- Eg) 
or, which is the same thing, 
Вес E 


(o3) - (2) 


Here м is integral, and » is essentially positive but not 
necessarily integral. 


875. The limiting form at which we have arrived at the 
end of the last article plays an extremely important part in 
the development of the general theory of Gamma functions. 
It will be very desirable for the student to pay considerable 
attention to it, and we propose therefore, in due course, to con- 
sider at some length the general behaviour of the function 


|J qd MES : 
ты ff 
(21) (z4- 2) (z4- 3) ... (# и)“ for different values of м and 


for different values of х, and the only restriction we shall 
place upon it at present will be that д is to be а positive 
integer, not necessarily large. 

Two theorenis, however, are required in dealing with such 
expressions as will arise, viz. 

(1) Wallis’ Theorem, which states that when m is a very 


pics d 2.4.6... 2n ; 
large positive integer, Е. Gat) and ~nr become in- 


finite in a ratio of equality, t.e. 


2.4.6... 2n yo 
Lty-© 173.5... (nI) nz | 
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(2) Stirlings Theorem, which states that when n is a very 
large positive integer 


1.2.3...» and /2n7.n".e-" 


become infinite in a ratio of equality, that is 


niet pa 
EE ul aeu 2T. 

The first of these appears in most treatises on Trigonometry, 
for instance, Hobson's Trigonometry, р. 331, Ex. 1, but scarcely 
appears to receive the prominence in the text-books that it 
deserves. The second, Stirling's Theorem, is less available for 
the student; hence these theorems are reproduced here for 
present use. 


876. Digression on Wallis' and Stirling's Theorems. 
ө? 


д°х? 


WALLIS Expressing sinÓ as выг) (1 1 to ©, 


and putting 9=5, we have 


:-6-36-36-9- 


1.8 3.5 5.7 (2n—1)(2n+1) | 


9$ sS gU Qn 
12.32. 5? ..c(2n —1) 


j Tm Ж ем (n) OE ee 


where e becomes indefinitely small when n becomes indefinitely 
large. у 
Hence, when n is TM we have 


2.4.6.. т eni 
RR E (2n4- 1) ultimately ; Ы 


and since n is very great, we have 


2. 4.. 
uc 1) ==, 
9.4... 2n. 


and 1:3... 2n—1) may be replaced by ~nr, these expressions 


being ultimately equal. This is Wallis’ Theorem. 


Its 
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877. STIRLING. Stirling’s Theorem states that for very large 
values of n, 1.2.3... n and /2nz e-"n" are ultimately equal. 


Write pln) for 1.2.3 ... п. 
Then $(2n) —1.2.3 ....... 2n 
and 2^6 (n)—2.4.6...2n. 
Hence Wallis’ Theorem, which may be written as 
25,49. 6?... (2n)* —- 
173/874 5. ez n anm 
gives tem mr. 
Let $Q). be called F(n). 
n” 9n 
Then 2?"[n^42n7 F (n) B—4/n (2n)?"4/ 4n F (2n), 
1.е. Е(2п) — [F(n)P. 
To solve this functional equation, write 2n for n. 
Then F(22n) = [F(2n) £ — [F (n) P. 
Similarly Е (2*n) — [F (n)? ete., 
and F(2?n) —[F(n)]*^, 


p being a positive integer. 
Now, putting 2?n—z, 


FG) Ее). 


Let p increase indefinitely and n decrease indefinitely in 
such way as to keep the product 2?n finite. Also let 


1 
Lt, ol F (n)]" 
be called Ж. 


Then F(z)=k*, which indicates the form of F to be expon- 
ential. We have to determine Л. 


Taking 1.2.3... nzg(n)—n*/2nmk^, 
change n to n+1. 
1.2.3... n. (n - 1) = (n4- 1) 2n 3- lrk, 


"uL ju. 
Hence, by division, А THR ee vanti, k 


n" P И 
ie. (1+2) (142) 


=e 
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in the limit when » is indefinitely large. Hence k—e-!, and 
therefore 1.2.3... and J/2n n"e-^ become infinite with m, 
in a ratio of equality, or, what is the same thing, 


е"! 
^5 d pers ам 27. 


This is Stirling's Theorem. The result will be considered 
further in a subsequent article (Art. 884). 

This particular form of proof was given by Dr. E. J. Routh 
in lectures at Cambridge (see also Dr. Glaisher on Stirling's 
Theorem in the Messenger of Mathematics). 


878. Illustrations of the Use of Stirling's Theorem. 
Stirling's Theorem is useful in such cases as involve factorials 
of large numbers. 


i. Thus the middle coefficient of the expansion of (1--z)^^ where x is 


! 
a positive integer, viz. es , is ultimately when л is very large, 
nm (9ny"e-? 92" 
пт mer oo aum 
This is the limiting form. It is of course infinite itself, but for large 
values of n a close approximation will be thus obtained. Thus, for in- 
stance, even T 4 case when x is not exceedingly large, in calculating 
! 
C= Go Y and uc from the logarithm tables the latter only exceeds 


the former by about 0°7 per cent. ; and in calculating C,)= От and 
9100 


Je the latter only exceeds the former by about 0:25 per cent.; and the 


error is diminishing as the magnitude of the numbers dealt with increases. 

Ultimately, for exceedingly large values of п, the middle coefficients of 
the successive expansions (1-- x)", (1+2) +2, etc., form what is nearly a 
G.P. with common ratio, 

92n42 / 92n 
grt on Улт?! 

аз is also directly obvious. 

2. The »^ number of Bernoulli, viz. Bj, , (see Diff. Cale., p. 502), being 
given by 


qu 4: T, 


2 (9n) ! 1 
Bonn = (2^ (1 gs gs 3n +.. c} 
we have, when n is large, 
2n. Эт (2n) e-?n 
(27 y^ 


=4т 244 e nnn, 


Bog 1 = 2 
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Similarly if E be the coefficient of z" in the expansion of sec + tan т, 
it is known dat 
gnis T 
Еа C ptt C pt peres, 


which embraces the cases of Bernoullian numbers and Eulerian numbers 


together, viz. 
К,„ = the n® Eulerian number, 


» 9?n(9?n — 1 
Kam ЕМЕ) В. 


(зее Diff. Calc., Art. 573, etc.), 
and we have when 1 is very large, 


. gne am 
Ka = ат тт nae” = sS e 


In this expansion, viz. 
UA 


sec v+ tan s= 1455 + T+ K+ — 
the ratio of the (n +1)" term to the n" is 
Šat 
Ки п’ 


а: — .n.-=— 
Té 


It appears that, since Lt z^ ==", the coefficients increase with great 


пі 


rapidity ultimately, and the series will be divergent for values of x 4 S 


3. In the series which gives rise to the Bernoullian numbers, viz. 


x a x 2% i: an 
goth, = ut g7l* В, 3- BEBE ...+(-1) ‘Вит! to? 
the ratio of the (п 4- 1)^ term to the n" is 


Leid uf ae 
By, (2n - 1) (2n) 
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and when n is large, 


= – Гі Agr 2d 6—2" 28-1 | а? 
= 4m ntie m+? (p1) (2n -1)2n 
iiaa {сс жаы 
Е se - D gs ras 
n 
UE A Е 
—-L^gbe£kaeU- v 
E 
~ dar?" 


The series is therefore divergent for values of x? + (27)?, and as 


Boy; (2n—1)2n n? 
ыч ERE ^^. ТИГУ z ultimately, 


the Bernoullian numbers ultimately increase with great rapidity. 


It will be noted that coth 5 becomes infinite if х have the unreal value 


2ur. When z is complex it is therefore necessary to limit expansion to 
the case for which the modulus of the complex is < 27.* 


879. А method of Calculation of the Numbers of Bernoulli and the 
Numbers of Euler is explained in the Differential Calculus, Art. 573. 
Both sets have been calculated for many coefficients of their respective 
series (see Proceedings of the British Association 1877), and probably far 
enough for all practical purposes for which they will ever be required. 
Several are quoted on pages 106 and 501 of the Differential Calculus. А 
few extra results are put upon record here for reference, for the con- 


venience of the reader. Also, as we are about to deal with such sums аз 


+ ... to © = Sp, which for even values of p are to be found from 


2(9n) ! 
Bua (rs Son; 


we tabulate a few of these results also. 
B,=}, B,— ds, B,— i5, B,— 3, В, =, By — d$, Вз=$, 
Вь= 389, By =y, By = PR; 
Е,=1, Е,=5, Е,=61, E,—1385, E,,—50521 ; 


T TS T? T19 
8-7, 51—00’ «7945» Ss= gag» S107 93555: 


The values of S, up to Sss reduced to decimals will be found tabulated 
later for purposes of evaluation of integrals to be discussed (Art. 957). 


880. For other methods of Calculation of Bernoullis Numbers etc., 
see Boole, Finite Differences, Chapter VI. 


881. We note that B, > B, > В, < В, < B,< etc., and the coefficient B, 
is the smallest of Bernoulli's Numbers, after which they rapidly increase. 


* See Bertrand, Calc. Dif., Art. 412, 
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882. The Value of II(). 
Consider next the case of Gauss’ II function for n=}. 


IL (5) =Le нй 


ue ow 
92.42.62... (Qu)? 


785 ZH 
3'àt 


esty 7 Д Өн. 5.0 on iler on а РИ 
Li.-«173 3.4... (39) Qu 1 T)" 
2 2% (u !)? 
wim Out ly” 
= 22 ит. ue. ei» ut 
— а 
= ІА, e : и 
& 1+5.) (+5 .) 24H 
( 2 ) 2u PAK 
$1.1 e. 
—eT sQU$3773) 
whence п(;) =“ ^ 


It will be remembered that for positive values of т, 
II(»n) = T'(n4 1); 


therefore r(3)= п(%)= Мк апа г(5)-5г(;); 


- г @=ут. 
which agrees with Art. 864. 


883. The Graph of y = 2"е-*. 
We shall next study the nature of the family of curves 
у=х"е® 

for various values of т. 

The subject of integration in the Gamma Function I'(n 4-1) 
viz. z"^e-*, has a maximum value when 

т2"—1е—2 —2"e7-—0, 4e. whenz=n (т> 0), 

and the maximum ordinate of the curve y=a"e~* for positive 
values of z is ne". | 

The graphs of the members of this family for n=0, т = 05 
ъ= 1, %=2 are shown in the accompanying figure for the first 
quadrant, which is all we require. 
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The case n=0, viz. у =e, is a logarithmic curve, and cuts 
the y-axis at a point y=1. It has no maximum ordinate 


Fig. 315. 


1 
The case т. = 0:5 has a maximum ordinate at x=}, viz. Же 
and then runs to the positive end of the z-axis asymptotically. 


The case n=1 has a maximum at z-— 1, viz. 2 


The case т = 2 has a maximum at х=2, viz. ah 
All the curves have the z-axis as an asymptote, and all go 
through the point (+ 3 ‚ Where they eross. 


For values of n between 0 and 1, the curves touch the y-axis 
at the origin. E 

Tl case n=1 touches the line y =g at the origin. 

The cases for n > 1 touch the z-axis at the origin. 

The several maxima, viz. n^e-^, diminish for various values 
of n from n=0 to n=1, and then increase again, all the crests 
the curves lying upon у —2*e-*, i.e. 


the least of the maximum ordinates being at х=1, and 
belonging to the curve у = ge, 
The area bounded by any of these curves y-—z"e-*, the 
g-axis and the ordinate at x=% , is 
[ew dz, 4e. l'(n4-1), 
0 
and increases without limit as « increases. 
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884. Extension of Stirling's Theorem. 
We have shown (Stirling’s Theorem) that when т is a large 
positive integer, 
1.2.3.... n 2 2n n^e-^, 
the meaning of the equality sign being that these quantities 
become infinite in а ratio of equality. 
We proceed to show that even when т is not integral, but 


still positive, T'(n4-1)2 2n ne, 
when n is indefinitely increased. i 
We have Г(п+1) = f z*e-* da. 


Let us transform this integral by putting 


жы. e 
z^e-*—qe-"e ? м 


which is legitimate, as n"e” has been shown to be the 
maximum value of z^e-7. 
Now, ав { ranges from — œ through zero to +o, 
х ranges from 0 through іо +o. 


Г(п+1) [* iede y 


Thus px PY X 
and we have to find 77. Let z—-n(14-7). 
Then (n EB nT) e ne-a = n"e "e LU ^ 


За 2 
5 (1+т)”е-"т=е 3" and 1юр(1+т)-т= -5. ...(2) 


Clearly т vanishes with $ and as Ё сап be expressed іп 
terms of т by expanding the logarithm, we can by the 
ordinary process of reversion of series expand 7 in terms of f. 

ций 
i prt Aeg ag 


Then, differentiating equation (2), 


Let T—À +.... 


d 
TA 0405 NEN TT PA e RRAN (3) 


whence, by substituting the sèries for r and equating 
coefficients, we can readily obtain the values of 4,, A,, A,, ete. 
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Now ама аг у eia == [| та 


mec dee 


© "а t 
=] e 2 [Art tag agi] 
" uu uu SUB. E 
and [me . dt = dw JT, 
by writing кё for 2 in the result of Art. 223, Ex. 4, 
| 2p+1 
г з) 
=, 
and | ttle—«'t* dt = 0, 
аз is obvious, for the negative elements of the summation 
cancel out the positive ones. 


Hence 
LL UMEN PR Lo) АГ AEG) + ote 
ый a 5. "Q 
2 Г. 2 
>. 2 = 1i 5 ( JN 


and it remains to obtain the numerical values of the coefficients. 
Substituting the series for 7 in the differential equation (3), 


(tyre art )* (4 As opi ) 
" (1+4 agp apt Ji 
whence м1 
ИА 
ba tao 
and generally 
1 CEST MET oats Hit + uu ia Dr 
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к т (n—1 —1)(n—2 
а Hi e arse ce i В шй 


+.. ‚+4 „Аз=т®А‚„.., 
А rh dy А 


14 * 9-3 


ie. (n+1)A,A en wo 
+...=пДА„ л, 
the series proceeding as far as the greatest binomial coefficient 
in (1+z)"*+1, and the last term of the, series being halved if n 
be odd. 
Thus A,=1, 
ЗА,А,=24;, 
4414,1-9.4,?—34,, 
54,44,4-104545—44;, 
64,4,4-15454 ,4- 10.4? —5A ,, 
74,44,4-21454,--35454 ,—6.4;, 
84,4,+284,4,+564,4,-+854А ,?— 7 A,, 


etc., 
giving A,=1, A,=3, A,=}, A47 — 5, Ag— ys, 
Ag— 185 4d — тў, A,— 11, A,=— 350150, ete 
Hence, finally, 
р. TI 
T(n+1)= J2Qn7 ne [Mz zal 


When a is indefinitely large, we therefore have 
(n+ 1)=J2n7 ne, 

which removes the limitation that n should be a positive 
integer, as supposed in Art. 877. Moreover, it will be noted that 
an expansion of eres is effected in powers of 1 , viz. 

C(n+1) _ 1+1 ELO B9 01. mni 1 
Jnn nhen E т 288 n? 51840 n® `` Ppl m? ^ "7" 
the law of formation of 4, +1 being as above stated. 


885. Ex. 1. In calculating 10! in this way, 
log ^/2т . 10. 1019 e-1? = 6:3561451 (Chambers! seven-figure logarithms) ; 
. 2 . 10. 101? 719 3598695 (the last figure doubtful). 
Carrying the series to four terms, viz. 
1+5 + 550 — rr dfggg = 100836537, 
we get 101— 3598695 x 1:00836537 = 3628799* etc. 
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The true value is 3628800, so there is only an error in the last figure 
in the approximation. 


Ex. 2. Calculate 100! Here 
log (100 !)=log {М Эт . 100. 1004 e-199(1 + зао — ---)} 
—157:9700036, 
indicating а number of 158 figures, beginning with 933262, viz. 
9:33262 x 10157, 


[The logarithms from 1 to 100 add up to 157:9700038, which is in 
agreement with this result, except for the seventh figure of logarithms.] 


886. Properties of Gauss’ П Function. 
We may now proceed to discuss the nature and properties 
of Gauss’ II function. 
Let us start again with a consideration of the expression 
А29, 
ее) CET 
where и із а positive integer, поё necessarily large, at present, 
and 2: is a fixed number, either real or unreal, positive or 
negative, integral or fractional, but finite. Call the expression 
II(z, и), and abbreviate it further into II(z) when in the limit 
u is œ, so that II(z) stands for II(z, oo ). 
Consider the graphs of 
| 1.5.8. 
У r3). GF)". 
for different values of ш: 
'ГВеге are и asymptotes parallel to the y-axis. 
y is positive from s=% to z——1, 
negative from z— —1 to z— — 2, 
positive from z— —2 to 2= — 3, 
and so on. 
And if и be >1, the z-axis is an asymptote at its negative 
extremity only ; 
also when 2=0, y=1 


when 2=1, y ii 
1.2,4* 
h = а 
when go V7 at DFE). 
etc. ; 
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and these ordinates approximate to 1, 1, 2!, 3!, ... as и in- 
creases, whilst at the same time the number of asymptotes 
increases. 

The cases of u=1, 2, 3 and 4 are shown in the accompanying 
figures, which are intended to exhibit graphically the general 
characteristics of the functions, but are not drawn to scale. 


The case »=1 gives E LAN a rectangular hyperbola, with 


z+1 
y=0, z— —1 for asymptotes. 
P 
о х 
Fig. 316. 
The case u=2 gives y=" ___ 9 
uere qvos YF FDF ` 


Fig. 317. 


- VAJ г, Те Arn v | 
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The case 4,—3 gives „д seri 


)(к-Е2) (ж-Е8)°` 


Fig. 318. 
S. ү it А 1.2.8.4 
The case и=4 gives V—irDGF3)SQRLIA) 
Fá 
= -2 =i о 1 2 3 
Fig. 319. 
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The lengths of the ordinates for various values of д and и 
are shown in the table: 


887. General Remarks. 

From these considerations it will appear that in these curves, 
viz. и=2, u —3, и=4, etc., 

(1) At z—0 all the ordinates are —1, and any two of the 
curves cross each other. 


(2) At х=}, 1, 2, 3, 4, ... the ordinates of the several curves 
form an increasing series, so that the curves as и increases 
are such that of any two the one with the greater u has the 
greater ordinate. 


(3) As z increases through zero the curves are all initially 
approaching the z-axis. The limiting case of the hyperbola 


juk continues to do so, the others all ultimately have 


2+1 
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ordinates >1, and therefore have minimum ordinates in the 
first quadrant. Moreover it may be shown that 
1.72 has a minimum ordinate between 1 and 2, 
и=3 T Е T 0:9 and 1, 
p= " d k 0:7 and 0:8, 


Аз м increases, the minimum ordinate begins to approach 
the y-axis, but does not do so without limit. For in the 
case u=% 16 lies somewhere between 0 and 1. 

(4) On the negative side of the y-axis at æ= — $ the succes- 
sive ordinates of the curves и=1, и=2, и=3, etc., form a 
diminishing set. 

(5) и=1 has one asymptote parallel to the y-axis, 

1.—2 has two asymptotes parallel to the y-axis, 
и=3 has three asymptotes parallel to the y-axis, 
ete. 
4. —1 is asymptotic to the x-axis at both ends. 
ш=2, и=3, и=4, etc., are only asymptotic to the z-axis 
at its negative end, and alternately from above and 
below the z-axis. | 

(6) Observe the behaviour between the several asymptotes. 

Between z— —1 and z— —2 the several ordinates at z-- — 3 
are all negative but numerically increasing, i.e. the more 
asymptotes there are the further do these branches recede from 
the a-axis. Similarly between the asymptotes z— —2 and 
‘z= —8, or any consecutive pair. 

Note also that for each given value of the branch between 
two consecutive asymptotes has a numerically greater ordinate 
midway between those asymptotes than is the ease for a branch 
between two consecutive asymptotes more remote from the 
y-axis. 

(7) The limiting case 

git. 
(z4-1)(z4-2) ... (z4- 14) 
becomes, when z is positive, the curve y —I'(z-- 1), as has been 
shown. 

The shape of this limiting form will be more carefully 
considered later in Art. 922. 


y 14, и“, viz. y —II(z) 
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But there is this difference between the functions 


AC ANS UIT: " 
= G4-1)(z4-2)...(z4- и)“ 
that though they coincide in value for all positive values of z, 
the former becomes infinite at the values z— —1, z— —2, 
x= — 3, etc., but has finite values for other negative values of 


z, whilst the definite integral is permanently infinite for all 
negative values of 2-1. 


Live and 7 
о 


888. That the factor form has finite values, when p becomes infinitely 
large, for negative values of х between the asymptotes may be made 
clear by taking а case. Take r= — $. 


Lu2;2... 


Then Lt,-« 


2.4.6... р 1 
1:1.3.5...( 2n — 3) n. 


а ln 92, 42.62... (2)? (901) 
1.2.3.4... (2p 302и 2) (8 1)(8и) i 


(J Qua pe ^Y (24-1) 


= – 49 ———® : 
м Au (25) e7 sje m 
_ 27E Qrp 2ң— z s. - Ns. 
9 mp и? 
5 23 
Similarly at х= -z the corresponding limit is I3 Мт, 

at z= —7 the corresponding limit ів e UM - 
PAS Оа R34 


and so on. 
These mid-ordinates, half way between the successive asymptotes, thus 
form a regular descending series 


S4 ug ЦА. $9. Бє 2 
-ivm 1:59 РМ Pes etc. 
889. It is worth noticing that II(z, и) may be written as 


| uy aT: 
П (2, DI CESITCESJTCES NN CEU D 


OG GENCE» y 
(1+2) (1-3) (14+§)---(1+2) uFi 
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(HOO HY у 


where E indicates that the product of all such fractions as 


follow " is to be taken from r—1 to "= и 
And in the limit, when и= оо, 


ta rs 


1+5 


ог, what is the same thing, when x із real and positive, 


0+) 


га +2)= P . 
^ (rez 
П 
890. Reduction of П (æ+ 1). 
Again, 
аа ава В RT 
FR) (@F 3) (23-4) ... Fp) @ Fe tI) 
=p Е П(х, и). 
Непсе 
1 
II(z4- 1, и) — (z--1)IL(, мх— т’ 


which is the law of connexion of the suecessive values of 
II(z, и) for unit differences in z. 
In the case when , is indefinitely increased, the factor 


ihe 
ЕБ 


becomes unity, and we are left with П(2--1)=(2+-1)П0 (a) 
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and changing х to 2—1, II(z)—zII(z—1) This is true for all 
finite values of x, positive or negative. 

In the case of values of z>0 we have II(z)=F(x+1), and 
therefore Г(х-+1)==Г(2), the formula already established for 
the Gamma function. 


891. The Case when wx is a Positive Integer. 
When z is a positive integer we may multiply the numerator 
and denominator of 
TUN 
| u= aM ub 
Me = ата) "era" "a 
obtaining in that case — II(z, и)= iw Be, 
and then removing и!, 
1:24 
П (2, 
© 107 G-ETY(Qa- F2) - ILS 
wm 


so that when м is indefinitely increased, x remaining finite, 
II(z) becomes z!, which is in accordance with the result 
l'(z4-1)—2z! of Art. 860. 


892. Comparison of the Gamma Function with Gauss' Function. 

It will now be clear, from Art. 887, that the two functions 
I(x) and l'(z4- 1) are identical for all real values of x greater 
than —1; but that II(z) is a more general function, embracing 
real or unreal values of z quite unrestricted as to sign. That 
II(z) becomes infinite for all negative integral values of т, but 
has finite values for negative fractional values of т, whilst P(x) 


defined as | e-" v*-! dv is infinite for all negative values of 2. 
0 


Graphically this means that the curves у=П(х—1) and 
y=T (x) absolutely coincide for all positive values of x, but do 
not do so for negative values of z. If we had restricted the 
definition of Gauss' function, viz. 
1.2.3. 
Ito Tos н) GET GER) GE) 

to real values of x greater than — 1, the identity of II(z) with 
Euler's Gamma function l'(z4- 1) would have been complete. 
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893. We have, from the definition, 
1.2.8...(и—1)и 
(2 п) —5y3—2)8—2) ... aE АШ 
1.2.8..(u—Dh ьа 
2(2+1)(2+ 2)... (2+0 — 1) " 
Hence multiplying them together, and assuming that 2 is 
not an integer, 


II(—z, 4) П(2—1, и) 
Mi 12. 22.32... (и — 1y? u 
2 (1?—a2)(2:—22)(82—2?)...((u — 13 —2?) их 
1 


А 
Tf, абам уы ав 
(1-2) (1-5) 41-27)" 


and II(z—1, и)= 


and when y increases without limit, Lt arh z being finite, 
and we have 
1 T 
II(—2)H(z—1)— 


It will be noticed that in proving this result no assumption 
has been made with regard to x except that it is not to be an 
integer, either positive or negative. For such values one or 
other M the II functions would be infinite, as also of course 


"would 
sin gm 


Taking positive values of 2 less than unity, and remembering 
that in that case П(2) = Г(2-+-1), we have 


l(1—2)l(z)—- 


sin gr’ 
as previousiy found. 

894. If we were to base the discussion of the properties of 
T(x) on this method of pu we could therefore Rod the 
value of the definite integral E me of Art. 870 to be = eret 
where 0 — z— 1, instead of investigating the integral frat and 


then deducing the result Г(1— 2) Г(х)= = 
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895. An Unreal Value of 2. 
We note also that if x be unreal and = vy, 
т . 
diC —g)1G9y— 1) x ту! 
but that Г, as defined in the Eulerian manner, loses its meaning. 
See, however, Art. 900 for an extension of the definition of Г. 
896. Both functions, viz. II(z) and T'(vx4-1) have been 
shown to satisfy the equation of differences 
Ue = (2-1) ,. 
Let us see from this point of view what can be ascertained 
as to the nature of the function uy. 
It has already been stated that this equation necessitates 
one form of the result to be 
м„= ÁAz(z—1)(2—2)... (т 1)ru,, 
where A is a constant or some arbitrary periodic function of x 
of unit periodicity, and u, is some initial value of u, to be 
chosen at pleasure. 
Following Laplace’s mode of procedure in such cases, assume 


as a trial solution, 
i w= tF (t) dt,* 


where the form of F(t) and the limits of integration are 
reserved for future choice. 
Then, since Ug, = (2+1) 2,, 


[ero dt=(c-+1) fero dt 
-[ro (2+1) 


рев) [enr dt, 


the integration being by parts, and the square brackets de- 
noting as usual that the term integrated is to be taken between 
the limits ultimately chosen. 

Hence the choice must be such as to satisfy the equation 


[e [F()+ РО] dt [РН 


* See Boole, Finite Differences, р. 257. 
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Let us then take F(t) so that F'(t)--F(f) —0, and the limits 
such that [F(t)¢*+1]=0. 

Our choice is now complete, and there is no further latitude. 

The first equation gives 79- —1, i.e. F(t)=Ce-', where С 
is an arbitrary constant as regards t. 

This determines the form of the function F in our trial 
solution. 

The limits must then be such as will satisfy the equation 


[Ce-tt*1— 0. 
Supposing 2-1 to be positive, this will be effected by taking 
t=0 апа t=o, for in each case Li -=0. 


Hence a solution of the equation for positive values of z+1 is 
== с| e~t dt 
0 


— OT (z-4-1). 


So &,—OT'(r4-1) is a solution, provided z+1 be positive 
where C is any arbitrary constant as regards t. 

To put the possible dependence upon:z in evidence call 
C, vz 


Then ` 1, v, l'(z--1), 
Wry = Vel (+ 2)= Vey (@+1)P (£+ 1), 
but Urn =(e+1)u,; 


г Us4177 Ug; 


whence it is clear that v, is either an absolute constant or 
some arbitrary periodic function of z whose periodicity is 
A 4- B cos? 2x Missa p. iue 
C+D sin? 2x bi aa байт ГИ? 
absolute constants, such functions returning to their original 
values when z is increased by unity. 

Thus w,-—f(z)T (z4-1) satisfies the difference equation con- 
sidered when f(x) is such a periodic function as described. 

It appears, therefore, that the equation w,,,=(%+1)u, 
is not co-equivalent with u,=I'(z+1), ie. Euler's Gamma 
function, or with w,=II(z), i.e. Gauss’ П function, but that 


unity, such as cos"27z or 
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these are particular forms of the solution, as has been pre- 
viously pointed out. 


897. Euler's Constant. 
The limiting value when n is made infinitely great of 


Leigh l-logn 
is finite, positive and less than unity. This limit plays an 
important part in our subsequent work. It is called Euler's 
constant and denoted by у. Its value has been computed to 
over 100 places of decimals (Proc. Royal Society, vol. xix. 
and vol. xx., p. 29). 
The first twenty figures are* 


у =0°577 215 664 901 532 860 60.... 
We shall presently show how it is to be computed. For the 


present it is sufficient to show that it is а positive proper 
fraction, and this admits of elementary proof. 


ы 
1 
ҮЗА ЗП И УВ mae 
аз tae pel a convergent series if 
r=1, 
ge 191 1 
=al g) (4—85) + 
— positive, since 7 2 1, for every bracket is positive ; 
M G41 og3)+(5 Hog 5)+-. EXC +1297) is positive ; 
b. 1.1 123 
- {ts 5+8 gt. .+- Elo O83'g'4'7 aj із positive; 
ГР 
eqptatagte +5 —log(n-+1) is positive ; 
and as log (n+1)> log n, 


i+ 2t 5+ te + log n is positive. 


*See Todhunter, Integral Calculus, p. 256; Serret, Calc. Intégral, p. 183; 
Legendre, Exercices, p. 295; De Morgan, D. and I. Calculus, p. 578. 
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Secondly, 
1 r—1 1 1 
zog = + log (1—7) 
1 1 "WM À 
== — 9,8 35-9. a convergent series i£r1; 


which, when n=, 


м > (; tiog oe "xa saa : "T m all convergent 


series, 
=a negative quantity. 


Therefore 

аа .+- * +log = + : са is a negative quantity, 
1.6. : + + deeds * logan i is a negative quantity, 
and .’. 34s s+; gt- ‚+- * tog ni is less than 1, 


and it has been ефе to be positive. 

Hence, making n increase indefinitely, у is a positive proper 
fraction. 

898. Closer Limits for y. 

Let tin =D 77 — log (n +1), va =$ r — log п (п> 1). 

1 1 

Then v, – w, = 106 (1 +1)- positive, if n be finite, and ultimately vanish- 

ing when n=%0, ie. и. =V» = у. 


з%—1 
rere negative; 


Now Un- Un = 1 +105 гы = positive; v, — Una -l +105 
tusrefore, as n increases, 4, increases and v, decreases towards the common 
limit у; and u,<y<v,, whilst n remains finite. 

Taking Bottomley's tables of Reciprocals and Napierian Logarithms, we 
readily find 

4,9069,  w,—'4014,...4,9—:5311, 4w44—:5532, 144,—'5610, etc. 
v,—1:0000, v,="8069, ... v,,—:6264, | v,,—:6020, V= 5938, etc. 
We thus have an approaching sot of inferior and superior limits for y, and 
note that it must lie between 0:56 and 0'60. It will be seen later that 

y =0°5772... (Art. 917). 


899. Except for negative integral values of z, II(z) is Finite 
whatever z may АН Real or Complex. 
Ни, Ua, Ug, ... и»... be any series of real positive quantities, 


each of which is М than unity, the infinite products I (1+4,), 


Ti (1— — и,) are convergent or divergent according as the infinite 


www.rcin.org.pl 


GAUSS’ FUNCTION IS FINITE. 89 


series Lu, is convergent or divergent (see Smith's Algebra, 
p. 423,* and Hobson's Trigonometry, p. 319), and if the quantities 
Му, Ug, ++. Ми... be complex quantities, the modulus of each being 


less than unity, the product TI (1--w,) converges if the series 
r=1 


x mod и, converges. (See Hobson’s Trigonometry, p. 320.) 
It can be shown that though the infinite product 


P(1+2), ie (1+2) (1+2) (1+2) (1+2)... to infinity, 


which occurs frequently in the present chapter, is obviously 
divergent, yet if we multiply the several factors by 


z 2 z 


e |, е?, e 3 etc., respectively,t 
we arrive at a product 


epe] 


which is absolutely convergent for all values of z positive or 
negative, real or и 
uw et 


is а series absolutely Аж if mod 2 <n for some finite 
value of n; whence 


" n) e intime 
= 1 Qm 
143 
n 
i) (1-2e* Lg mt) 
.е. ^ 
-1-£. s), вау, 


where e, is а series absolutely convergent which for finite 
values of z ultimately vanishes when x is infinitely large ; 


Р(1 tZe*-Ph-5 UAI Le) 


* Also see Arndt, Grunert, xxi. 78. 
TWeierstrass, Abhandlungen Acad. of Berlin, 1876. See also Hobson, Trigo- 
nometry, р. 327. 
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Suppose Ё the greatest of the moduli of 1+, for all values 
of z within a range for which the greatest modulus of z 


does not exceed a given finite quantity, then p) mod DP 
is an absolutely convergent series, and therefore also 


© 2 
D allte) is an absolutely convergent series, and since 
1 


Pü Чи,) is absolutely convergent when È mod v, is convergent, 
1 


Parae" 


is an absolutely convergent product, as is also 
© 2 £ 
Ри) 


Now Gauss' II function being defined аз 


142,2. z 
пе, ое) FH) 
и? 


can be written | —Lt,-« 


(0633-02 


where у is Euler's constant, which shows that for all values 
of 2, real or complex, positive or negative, excepting negative 
integral values, 
e Y 
IGI finite Fonction of 7' 


and is therefore finite. 


900. Extension of Meaning of Г (z). 
So far it has been convenient to adhere to the Legendrian 
definition of the symbol T (æ), viz. 


Г(2)= -[* e-"v*-1 dv, 
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and to regard « in this Eulerian integral as representing a real 
variable. It has been shown to be identical with Gauss 
II function, II (2—1), for all real positive values of a Having 
drawn attention to the difference of behaviour of the function 
defined as an integral and the factor-function of Gauss for 
negative values of z, it is scarcely worth while observing the 
distinetion further, and we propose to extend the use of the 
symbol I'(z) to negative and unreal values of z, which means 
that, when 2 is negative or unreal, Г is defined by 
n "T 1:9 лди 
РПО В, Тен ... (6 2) к rss 


and that when z is positive it is defined either in this way 


2 
› 


ог аз | е-°1 dv, and therefore we shall in general regard II (2) 
0 
as identical with T'(z4-1) or zT (2) for all values of 2. 


901. Thus а meaning will be given to such an expression 
as I'(a+./—1)), viz. 
9 pated 
Muna at) а тат) 
EE GE (HP). е) 
e- v) 
~ a finite function of (a+b) 


(Art. 899). 


902. Ex. 1. The modulus of I'(3--:a) is МГ +:а) Г({ -«a) 


=N {T (14a) га -peaj-A (Art. 895) 


sin ($ +:а) т 


T 
cosh a ' 


Ex. 2. If 1, a, o?,... a" be the л“ roots of 1 (n odd), we have 
(1+4) (1 +а2) (1 +a)... (1--a-12) — 1-2, 
and 1-+а-+а?+...На"1=0. 
r=n-1 
Hence П(х)П(ах) П (оёх)... (а) = Р II (ax), say, 
r=0 


r=n—-1 Za" 


pn 4, trem TEE SWR PET a 
(1-55) m) 0 
1 
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Б (1 +5) (1 +21) (1 t5) 5. to ut т т eS) ose Elaa) 
За bn da sul Сор Мы x 
"Р (Iz) P Г(1+ az) Em 

thus (1+ Se Fern era d ollis 


where 1, a, a?, ... are the n™ roots of unity. 


903. Gauss' Theorem. 
This theorem is a generalization of that of Art. 872, and 
includes it. It states that for any value of z 


nT (2) 1(2—>) 2—5) mc) 


"1 nd 
Tr) SML г а 
ог, what is the same thing, аз will be seen, 
1 2 —1 
тегә) T (2) rz)... Г(2+" n-1 
= (2«) ? nb. 


T (nz) 

Let the left-hand member of the first equality be called 

$(z) Then, first, we shall show that ¢(z) is independent of 2. 
By definition, 


muy Le) | Mrs dod А" 
п, 2 fes re 2) (ets) 


he 


1 ии —r) di 7. 


ita +nz— E 


n-1 
i nZynk&,nZz, 2 (п 
OR aut Cae GU 5 al 


1 
Е пп) (2—7). (2— ^; 
where D is the product of the factors 


n--nz, 2n--nz, SURE, | .—- un - nz, 
n+nz—1, 2n+nz—1, 3n+nz—1, ...... un A- nz — 1 
n-J-nz-— 9, 2Qn+nz—2, 9n--nz—2, ...... un -nz — 2, 


n--nz—(n—1), п-т (n —1), 3n--nz — (n— 1), 


t.e. 


[(nz+1)(nz+2)... 
—(nz4- 1) (nz4- 2)... 


(nz4-n)][(nz 4- n 4- 1) ... 
(nz+ ип). 


.. un+nz—(n—1) 


(nz4-2n)]...[... (nz -F ит)] 
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Hence 
»non(z-1).n(-523) - n mmt Чыр. 
Again, writing mu for u in Gauss’ expression for П(и2), 
Io) Ii ce Sehr cy 

soa сы З шу 


EIC. 


(nu!) 


Hence 


paw N" ug 


from which the z has disappeared. 

Hence, $(z) is independent of z. It remains to find its 
value. То do this we may either obtain the limit of the right- 
hand side directly, or avoid this by comparison with a known 
case, for а particular value of z, which will be a legitimate 
process, inasmuch as its value, not containing 2 at all, is an 
absolute numerical constant containing n. 

Adopting the direct method and employing Stirling's result, 


n-i ры, n 
ео пни" Dum ute Н 
J nee pe T 


=] 
nup (2m) 7 зиме _ (2x) 


pind УЕ n? 


| “| 


= It 


Hence, finally, 
»"Iro)n(s-1)n(s—)...m(4— 5.) Ont 
Teen 


904. If we adopt the plan of comparison with a known 
case, take the case of a real value of z, viz. 2—0. 
Then, remembering that П(2):=Г(1 ў 


o(e)=9(0) = rayr(1-5)r(1—2)...r (9—7) тга); 


or, reversing the order, 


-r(jr(9r(3. т. by Art. 873. 
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Writing II(z)—T' (z4- 1), etc., we have 
пг 1 8. 
"Ге рен”, гея) („уг 
l'(nz4-1) ni 
ie. reversing the order of the factors in the numerator, with 


the exception of [(z+1), and writing ['(z+1)=zI(z) and 
I'(nz4- 1)-2nzT'(n2), 


nare) (2+7)r (2 +2)...P(2+"—) NOS 


— 


nz (nz) nè 
nener (а т)г(а+ =)... г(а+" =") na 
Ry eer Я п 


which may be written as 
г(@)г(@+5)г(@+)...г(а+"—)—г(ш)(зх) F (hn. 


905. Cases of Gauss' Theorem. 


Putting 2=1 we have the result of Art. 873, viz. 


(rr). түе) вн 


Particular cases are 
n=2, Г@Г(т+ +5) =P(22). (27)523-", 
ie. r()r(z4-D)-,z reo, 
ай 


Гог 2, 
Керген 


n=3 gives го) г (2+1) (2+5 = 3s- sagr (32), ete. 


ie. putting 


906. The case n=2 may be deduced directly from 
4-1 9+1 
rrt ) 


| sin? @ cost 0d0— - 
| к) 
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For putting g=p, we have 


А PEDI 
[s 70 po „д ЛУ, 
А n? 0 cos = ГО ; 

. p+1)\\? 
. P 
у» ЭГЕ. 


and writing 20 — 9, 
ү; 311220 40 =f sin? ф dg- |. віп? ф dp 
T)re 
(E) 
Безу rege 
TT) (27) 


ie. WT (5) Г (PT гр 1). 


907. An interesting proof of this result is due to M. Serret, 
(Calc. Intég., p. 174). 


1 
Since B(p, o-| 2 (1—2)‘—14х we have 
0 


B(p, p)=| &—25714:— [i 8a. 


And since the integrand assumes equal values, whether 
we put «= 4-h or 1 — Л, its values are symmetric about 2= 4. 
Hence 


" "s 
B(p, p)=2 1—(@%—*Р-\&. Writing — , 


Вр р) | ar 0-27 (7 Е) 


1 f! 1 
== ail 7 (1—2)»31dz— 98i BG, р), 


N о or 2:3 (р)Г(р-+)=МяГ(2р) 
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or writing 2p—q-4- 1, 
9+1 а+2\_ - 
xr (LE) r(£25)-/zrq--1. 
908. Another form of the general theorem is (wri ting = for г) 


Or Er ERr -roen a 


N n 


ie ГГ (2)... ге") гг + 2) (2m) т 


n n 


2+1 
909. То ргоуе | log T (x) dz—2 log 2 — х- $108 27. 


zr 


Taking Gauss' Theorem for a real variable z, 


2 2 a-i 
P (9)T(z-- )T (2+2)... r (a) (ne) (2r) T ni-n, 
we have, upon taking logarithms, 


1 sci 
zlog Í Г(па)(®т) * nin 
] 1 ®—1 
=: [log T (2) +logr(2+2)+...+logr (s) 
=», : log r(z--7) from r=0 to r=n—1, 
1 
= logT (2+ y)dy, when n is indefinitely increased, 
0 


+1 
ый, logT (v)dv, if v be put for z--y. 
» z 

Thus, by Art. 884, 


E: ы Z minn 
n © L = (2) nie 


dx NI о — NX n—1 
| { log '(v)dv— Lt, 1 g Гу тат (nzy"e-" | 


=} log 27+2 log r—z-—log a* e-* (Эт). 


910. This expresses the area bounded by the z-axis, the 
curve y=logT (x), and two ordinates at unit distance. 
Changing z to +1, and adding to the former, 


z42 с 
| е log l'(z)dz—log (a* (x --1)* 6727 (т), 
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and so on, and more generally, 
3 
| log ['(x) dx 
z 


—log[s*(s-L-1y* (ed 259... (2n — eint 
where n is a positive integer. 


iHe 


? (my, 


-m"nr- 


911. Expressions for the Differential Coefficients of the Function 
y(x), log l'(2--1), and Expansion of log Г(ж-1). 

f d T(x) 

Let us write \/ (2) for $ log Г (x), t.e. T) 

Then taking the logarithmic differential of Gauss’ Theorem, 


(nz) и” ГО) P (2+ =)... Г (e+ ey | Qr) = nt 


ny (nz) =n log пануе)... (2+9=1) , 
апа differentiating again, 
, ; , 1 ; —1 
n*y/(nz) (а) + (+ 2) (a). 
Непсе 


пу (nz) = У) = (2+7), from r=0 to r—n— 1, 


ie La, mna) = [oy v err. 


=y (£+ 1) уве ове). log Г (2) 
= og р mazo в, 
ie. Lt, 4 (nz)y/(nz)—1; or writing v for nz, У in the 
limit when ~v is infinite, and therefore y(v) ultimately 
vanishes. 

That is = log T (x) vanishes when z is indefinitely increased. 


da? 
R l'(z4-24-1) 
Now Геи) 


Hence, Lie the logarithmic differential, 


1 1 
ЫЛЫН eg 21 z42 ^C ниен) 
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and differentiating again, 


nO NET 1 1 1 l 
Yu mou cape i 
and it has just been proved that y'(x+n+1) ultimately 


vanishes when n has been pron increased. 


n Fog Гав 


The series (1) is obviously convergent for all values of х > 0 
becoming infinite at z—0. 
Integrating this equation between limits 1 and a, we have 


ve-vo-|- Hz 
-(i-2* G- ze) G-z:3)* me rine (9) 


which is a convergent series; for the test expression, viz. 
n(x+2n) —9 
mathet 
and is greater than unity. (Зее Smith's Algebra, Art. 342.) 
Again, we have seen that 
1 n—1 
ny (nz) =n log n-- y (x) J- V (s HD Hye (2+ - ) у 
and putting z— 1, 


y (n) —log n4- У) ly (1+2), from r=0 to r=n—1. 


- "HR iSt Bed ra 


Lt. ( 1 Mes) pe 


Hence when т increases indefinitely, 


І4,-„ [ү (n) —log =] уа) dz 


=[овг(+ | =т= jog 10. 
diia. Ld „(КЕ - орт) =0. Lee (3) 


Putting = in equation (2), 


| Wo )—V (1) 1-5 3-- to 00 , 
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i.e. by equation (3), 
— ии, (++ +2 Лор п) 
— Euler's Constant у, 


; d 
1.6. y. (1), or I5 log T) oim na RENI) 
Hence, by equation (2), ti 


т o =)+(5— ait to oo 


5241 ———7 + ... t0 %0,......(5) 
which may also be written as 
£ log T (x+ 1) = Lt,- орт у I. 
ах E 214-1 2-2 z4-n 
Again, differentiating equation (1) n —2 times, we have 


а (—1)"(и— n| 4+ + кыз+ < | 
(6) 


$6: NI (Lor (тет. 1 =(—1)"(n—1)!8,, 


where S, deb as vg e 


which is convergent if > 1; or, what is the same thing, 
{ок r+} _, ад, (T) 


jp fiog r(æ+1)) —logT(1)—0; 
=0 
we thus have С 


{log Г(=+ у =0; FAL T'(z4- n} - i; 


and [елок — —(—1)^(n— 1)!S,, where n is 4 2. 


Maclaurin's Theorem "98 gives 
а? 23 4 n 
log T'(z--1)— —y24-8,7 —$,5 S, — (1S, Tes 
a result otherwise established in a subsequent article, and which 
will be thrown into a more convergent form, by the addition 
of other known series, for working purposes. This series is 
convergent if z be numerically < 1. 
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912. Collecting for convenience other useful results of the 
above article, we have ` 


(а) Lt, s 72198 Г(@)= 0 and Lt ишге, апа їп 


d 
any case 75 log T (x)= Е to oo 
and is positive. 
T'(n) _ "3 
(b) Tin) =logn when nis infinitely large. 
d 
(c) tog Г(=+ DEF Ferg and .. teg re --» 


(d) ор Г(®-+Е1)= —+(— zu)*(G-z:3)* ——M— T 
(e) Since Z log T(x) is continuously positive for all positive 
d , Е | : 

values of т, 25198 Г (2) is an increasing function as x 
increases from 0 to oc, starting from the value —oo 
at x=0; or, putting this geometrically, the tangent 
to the graph of y=log T (x) is continuously rotating 
in а counter-clock wise direction as ж passes from zero 
to infinity; and the curve is always convex to the 
foot of the ordinate. 


913. The student may note the following particular values 
2 
of Z log T(z), te у), viz. taking л?=9'8696044011, 


y (0) Q9, 
y'(5)— i eror над. J- тт — 19948022, 
9 
т? 
y'Q)- + s =1:6449341, 
у 1 ae RENI R 
A R = 9348029, 
9 1 Es i 
у ( у ааа +. =! = 6449241, 
‘ 1 TAI BV ui 
Аа ce i ij = 4903578, 
bd g* d Lp 
y'(3)— = tgtpgnt-78-pgp-g$7g-l 25 = 3949341, 
ete. 
У (<) =0, 
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which indicate how Z log T(x) is decreasing as x increases, 
but always remaining positive. 


2 
914. Since log Г(2+1)= —y2+8,5—S,5+5,G—--. we 
may write I'(z--1)as 


zo а% 


l'(z4-1)—e* ^O" er б 
=(1 ну aa oo Be c1 


: zi 
=1- ye (uS) 5 0 -3y$ 38). + 


which expands Г(х--1) as far as cubes of 2, and which might 
be useful for very small values of x, but the presence of 
powers of y renders calculation troublesome, and less incon- 
venient methods of calculation will be given later. 


915. It is noticeable, too, that 
log Met) урна grip _.. 
and that the several differential coefficients of this expression 
are therefore free from Euler’s Constant y, viz. 
d^ log l'(z4- 1) 
ах” x 


ых S, n42 (n T- 1)! » (n+2)! 
=(—1) 1{ 8 Set pj Said ED ма НИХ T En 


EN n— m 1 n+l Rn+2 п+1п+2 8,43 - .. 
Pic (a El быз 1.9 9+3“ ep 


And, similarly, if m be any positive integer, 
d^ dy Чья. WR 
Jat” log P(a+ 1)= (7-) [-ve Xie ] 
© S, 
= — (m4 1), yenit E тна, 


where (m--1), denotes (m+1)(m)(m—1)... to r factors, if 
n <т- 1, and is free from y if & 2 m-4-1; also that 


per r +1)| = — (m4-1)! 
нш" ок Ге) | _ =—(т+1)!у 
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916. Expansion of log Г(1--х) deduced from the II Function. 
The series E А a 
log Г(1+2) = —– ух+8, 5—8: 3+8 4— ud 


may be arrived at at once by taking the logarithm of the 


Gauss formula in the form 
а 


T(e-42)-Lt,.,——————. 
| np 02 
viz. 
log P(1+2)=alog u—log (1-- 1) —log(14-2)—tog(1--5)—...; 
and expanding the logarithms, supposing —1 < z < 1, 


2 
log '(14-z)— Lt | x (log 1-8) 8,7 8,74 Да 


ЧИ: ИИТ! 
where 8,= На +... 


and Lt(S,—log и) = Euler's Constant y, and the series S, (r > 1) 
are all convergent. 


Hence, 
2 n 
log (14-2) — —Yyet 8 8,2-8,2 ы +(—1)"8,——+ UN 
(—1<2<1). (1) 
Now, the even terms may be removed by the addition of 
jlog Si. 
sin zr 


For sin Sr (1-2) (1-2) (1S) ad inf. 


and taking logarithms and expanding, 


de zT 2° x 
O= 3 log Sg a RTT (2) 
Adding to equation (1), 
x x 25 
log Г(1--2)= 3 log meg yt Se 3 Ss BT eB) 


The coefficients S,, S,, ... all begin with a unit. This may 
be removed and the series reduced to a much more convergent 
form by the addition of the series for tanh^!z to each side, 


Vız. 1 +2 


tanh-!z—4 log tos T ET ^e 
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And we then жеке. 


log I'(124-2)— th “——tanhz-+(1—y)« 


sin 2 
 -&-nZ-&- n2. iino nj 


The values of y, Sz, Sg, ... 9 are all calculated, and the 
tabulated results are given in Art. 957. Euler calculated 
S, to Sı Legendre* gave the values S, to S, to sixteen 
decimal places. The list in Art. 957 is taken from Legendre's 
list as given by De Morgan, Diff. Calc., p. 554. The series (4) 
converges rapidly and is used for the calculation of the values 
of 100 Г(х). Legendre gives a table of values of LT (x), t.e. 
10-106 T(x), from LT'(1:000) to LT'(2:000) to seven decimal 
places, in his Exercices du Calcul Intégral, pages 301 to 306. 
A table is also given by Bertrand, Cale. Int., p. 285. 


917. Calculation of Euler's Constant y. 

These series may be used for the calculation of Euler's 
Constant y by taking a value of z, for which T'(z) is otherwise 
known, viz. z—1, for which T'(z)— 7. 

Equation (1) gives 


1 4 Сау ! 
у= mi Г(2+1)+5,5—9:5 8,1 — 5 
and putting х=}, 


Á Let Zar adar | 
y -log, 38; à 3959 44 98 - oree өөөөөөөө (5) 
Equation (3) gives, by changing the je of z, 
log Г(1— 2) = } log -- T = yr 8, 3 xg e B z+.. 
and putting х=$ in s 
ar oem T C 
y log 2— Si BSsgi Өт oo ve ee ee (8) 


which is more rapidly convergent than the former. 
Formula М gives 


3 ЕН 


log gy -logT— Plogg4 5 Заа 3 9 B 95 
И де 8—11 5—11 5,—11 7 
4.€. у= log. 5 ——3— 5? SE M 34 7 NN eee € M uv 


* Traité des fonctions elliptiques, Legendre. 
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This is the best of the three series to employ to find у. 

And with the aid of the tables of values of S, the calcula- 
tion to seven places, which is all that is likely to be wanted 
for ordinary purposes, may be readily performed. 

The value of у is 

y='57721 56649-01532 8606... , 
and ‚1—у='42278 43350 98467 1394.... 
The value of log,10 is of course required. It is 
log, 10=2°30258 50929 94045 68401 79914... , 
and the modulus log,,e='43429 44819 .... 


918. The numerical calculation of values of log I'(14-2), and 
therefore of T (æ) itself, will now present no difficulty. With 


the values of saul 4 s s: etc., inserted, the working formula 


stands* as 
log T-c2)-j log T og ot 42978432 
‘aaa тт 2 
— 0673523023 
—007385525 
—001192727 
— 00022312? 
— ete., 
and is rapidly convergent for the small values of z less than 
z=}, 2? being 1024. Hence the last term :00022312? in 
the case x=} becomes ‘8000004, whilst for z—1, which is the 
largest value of z for which it will be necessary to use the 
series (see Art. 921), the error in omitting all the remaining 
terms of the series will not affect the seventh decimal place. 
Hence we have here all that is necessary for the construction 
of seven-figure tables for log T (x). 


919. It is worth noting that the addition of log (14-2) 
and log (1—2) respectively to Г(1--2) and T'(1—2), viz. 


2 4 
log I'(12-2) = — yz4- S, 2—58, pe Su. 


2 4 
апа log Г(1—2) = ys 8,7 +8, 5-8, «Te 


* Bertrand, Calc. Intégral, p. 250. 
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give logI'(14-z)— —log (1+2) - (1— y) 
: 2% zs qt 
+(S—1) 5-——(5—-1)3 +(S,-l) gT 
and 106 Г(1—2) = —log (1—2)—(1— y)x 
gi a? Р аА 
+(%—1)®+(5—1)--+(5,—1)--+...; 


whence 
4 log Pu - — tanh-!z4- (1— y) z— (S,— 1) Е 
a? 
p CRINE 
But flog Г(1-+2) Г(1—2)= blog 22 2» i.e. adding, 


log T (1+2) )=} log 57 Z-—tanh-tz+(1—y) 2 


© т 
z^» (Sensi — Dau ru 
the same series as before, which шау be written 
waz l—c а gent? 
log P(1-+2)=$log (— T 155) ty) ey ne xm 2n41 7-7 UTi 
| ; 1—2 —1 1 
and putting z—1, since Lt,-, pet opt ods to mag 
Sensi — 1 


1— утаа ; 
and putting х=, since T'(3) nm 


M eem 
These series are given both by Serre and Bertrand for the 
calculation of (1+) and y. 
The formulae 


1 TE 1 1 
log P(t +2)=5 log BER у2— 5832 — = S52°— ©, 
1 1 
log Г(1— —1)=5 log ds Буга ls A RE 5552... 
and ye I Eh PE. . 
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were given by Legendre (Exercices, p. 299). But the addition of 
the series for tanh-!z adds to the P of the convergence. 


920. Since T(m) '(1—m)— , we have, on putting 


sin MmT 


ET for m, 
1-х 1—zY _ T ron ; 
FF) 0) 
sin T cos 
2 2 
But га Je Г (Art. 905). 
Г г(1+=) 
Hence, writing 2 in place of z, 
г(2) 9-60...) 
EE 
2 4 
From equations (1) and (ii), eliminating rh), we have 
z 
ШГ r(2) Ы 
Г (2) = —— —— 06060) 
21-2 сов” D (525) 
2 2 
921. By means of the four formulae 
Г(2) -(z— 1) l'(z— 1), ...(1); Г(а)Г(1—2)= 7, HER AA (2); 
| ЛА 
a /— Г(22) JT 
Г (2) =21-2,/т —1——,.(3); T= ———— ‚ (4); 
| Г(+2) "Y сов 7 Е Tí à a 


it may be shown that T (x) can be calculated for all values of 
x when those between Г(1) and l'(3) have been calculated. 
(a) For 1< 2 < ә, reduce by continued application of 
formula (1) to a case 0 < y < 1. 
(b) For $ < 2 < 1, reduce by formula (2) to a case 0 < y < $. 
(с) For + < < < $, reduce by formula " toa case } <y < 3. 
For if 25» 5 3^6 and t e 


and if 2« 3X3 and > 


8 
m 


3” 
i 
ó 


WWW.rcin.org.pl 


GRAPH OF T(x). 107 


(d) If < x < 1, the case needs no reduction. 

(e) If 0 <= < 1, use formula (3). This involves Г (34-2), 
and $ х lies between 4 and 2, and therefore falls under case 
(с), and an application of formula (4) reduces I'(z4-3) to cases 
in which the arguments lie as before, viz. «y < 4. 

In Г(25), which occurs in the numerator of formula (3), if 
0<2< 1, we have 0< 22 < 1, and И 2z ^ 1, no further 
reduction is necessary. 

But if 0 < z < 15, we have 

0<22<1 and 0<4&<1. 
We then use formula (3) with 2x written for =, 
: 2а ү 2) 
е. Г(2хж)=/т2 TEF) 
Similarly if 0 < £ 34, use 
dee 1—8 Г (82) 
Г (42) =/т2 Tatia 


and so on. 

Hence it follows that the use of series will be only 
necessary in the case of I (s), where z lies from 4 to 3, and 
that when this group is calculated by the series, all others 
follow by the above rules. 


922. Graph of y 2 (x) -| e7*2*- dz. 
0 


Regarded as defined by the integral, it is plain that so long 
as z is real and positive Г (x) is a positive function, and that it 
becomes infinite if z—0, as may also be seen from the fact 
that Г@)= Г (2-1), and therefore r()- 00 36. 

We have seen that 
d? 
dz 25 log l'(z) = tui (24-1) = gt. 
and therefore is infinite E z —0, but for all values of x 
from 0 to со it remains positive and finite. Hence 
d eh Pte) 

di log Г(2), ùe. l'(z)' 
is an increasing function of 2, and its value at z—0 is 
obviously — оо, for 


2 tog T (@)= -y -=)+(- zn) ... (Art. 911). 


— sd .—— 
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.Also, when z is +0, 


рун А. to о = | о. 
H Г(2) . 
ence q (2) increases from — æ through zero to +% аз x 


increases from 0 to oo and as T (x) remains positive throughout, 
T(x) changes from negative to positive once, and once only, 
as ж increases from 0 to х. 

Therefore T (x) has one, and only one, stationary value, and 
that is a minimum, and I(x) decreases from o» when z—0 to 
Г(1) =1 when 2=1, and since Г(2)=1 and T'(1)— 1, the ordi- 
nates at z—1 and z—2 are equal, and the minimum lies 
somewhere between z—1 and 2-2, and is numerically less 
than unity. From 2=2 to x= the value of I(x) is con- 
tinually increasing. 

The curve then 

(a) lies entirely on the upper side of the z-axis; 
(b) it is asymptotic to the y-axis; 

(c) it has a minimum between z —1 and z—2; 

(d) it recedes from the z-axis from v 22 to == о. 


The equation to find the exact "n of the minimum 


ordinate is 2T) o, or writing z —1-Ff, S Г(1-+Н)= 0. 
dlog (1+4) I'(1-t) 
М ЖОЛГЫ 
fiance T Г(1-+Н)= г(1+9| — p y)H-(,—1)t 


TE (S,— 1)?4- . zt 
and ¢ is to be found by trial from 
1 z 
T= 0422784 ... + (S,—1)t— (8,1) 
and substituting for S, and S, their values in decimals to a few 
places, an approximate value for ¢ may be obtained, and by 


the usual approximation methods the result may be found as 
nearly as desired. Serret gives the result to seven places, viz. 


t=0°4616321... 
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i.e. the abscissa of the minimum ordinate is 
z—1-4-:—1:4616321..., 
and the value of the corresponding ordinate is found to be 
y —0:8856032....* 


In the tables for LT (x), i.e. 10+logT (x), we find in the 
vicinity of the minimum 


z L Y (o) $ LI (x) 
145 9:9472677 1:463 9:9472396 
1:46 9:9472397 L4: 9:9472589 
1461 9:9472393 148 99478079 


1:402 9:9472392 


RANE NUTUS 
е ЕД е S 

ы ДЕ ы р е Г 
Е Ета Ы Е 


и 
Е 
a 
a 
= 
* 
: 


о2о 1111 
CESNN эче 
а ЖШ Ж ЕЙ ЕТА 
ЕЕ 
Е Е: 1328] 
obleko] | | | | 


So we see from the tables that the minimum ordinate is in 
the vicinity of 1:462, and the value of the corresponding 
* Bertrand gives 0:8556032, page 283, and again page 284, line 3, and the 
result is given elsewhere. This is evidently an error. The result is given 


соггес у in Serret. Cale. Intéa.. n. 186. Suiv ao prre 
ггес у in Serret. Calc. Intéa.. n. 186. d 
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logarithm, 1:9472392, indicates an ordinate 0:885603 approxi- 
mately. The minimum ordinate is reached, therefore, a little 
earlier in the march of z from 1 to 2 than the half-way 15, 
which might have been expected from the very rapid fall of 
value in T(x) between Г(0)=о and Г(1)=1 and the much 
slower rise on passing х=2, Г(2)=1, Г(3)=2, Г(4)=6, 
T'(5) = 24, etc. 


Г(24-1) JIar ate 
z 


For large values of z, approximates to EC 


and the graph of y=T (x) to the curve y= T (y. 


We have now seen to what shape the several curves in 
the graphs in Art. 886 are gradually tending, and com- 
parison should be made between the figures given there and 
the graph of the limiting form у= Г(х) in Fig. 320 of this 
article. 


923. It will be noted that since T(x) is decreasing from 
z=0 to z—1:4616321... and increasing from z—1:4616321... 
to =o much more slowly, the differences are negative for 
the first part of the march of I (x) and positive for the second. 
Similarly for the differences in the tables which give log T (x) 
ог LT (x). The tabulation is only effected from z—1 to х=2, 
for by virtue of the reduction formula I'(z4-1)—zT'(z) this 
is all that is necessary. In using the tables care should be 
observed with regard to the change of sign of the differences, 
and those who wish to make close calculations should observe 
the remarks made by Bertrand, Calc. Intég., p. 284, with 
regard to the behaviour of the differences both of the first 
and second orders. 


924. The rule of interpolation commonly used is 


2 


us m ty bx Au О дц... 


(Boole, Finite Differences, Art. 2), 


rather than the ordinary rule of proportional parts, which 
stops at the second term. 
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925. Expressions for 


d 
ia log T (x), 
as definite integrals. 


d? d” 
T — log Г (2), Lis etc., 


2 
The expressions for gi log T (x), 2 log Г(5), etc., viz. 


d 1 1 
ВТР Ltn- - {logn- (tz xi aho 


d? В 
dà log Г) тур lh ques B um (2) 
T log DI'(z) -(— 1)" r(»)| + toi eim . to «| (3) 


can readily be converted into definite integrals by aid of the 


lt 
results [ie e-&g-1dg— s ho Neck саім; (а) 


© —t. p—kz 
ind | CEU dbi сойлой odd. ner (b) 


0 
(a) has been proved in Art. 864. 
(b) can be established thus: 


=k? 1 
ва | =>. 
J. wide k 


Integrating with regard to k between limits 1 and k, 


оо —kz k © 7—2: p—kz 
771284) Ee ] dz=f ate: 


To convert 


d E H 1 1 1 
25198 T2) = Lt, {log НЫ vasi 5 рая) 


the right side may be written, by aid of (а) and (b), 
=һ-| | [Ar ree e aei iis ‚—е-ва+"-) ав] 
0 


B 
=Le | ee A ) ag | 
ам 15) == et aja 
-| Cg- нь (А) 


for the second integral disappears when n is made infinite. 
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926. With regard to /у =| e Tn s estt ri it may be 


desirable to make a closer investigation, for though for all 
values of В between e and infinity where e is a given small finite 
quantity the factor е" destroys the integrand when n is made 
infinite, there may be some doubt as to the behaviour of the 
expression in the immediate proximity of the lower limit. 

We note that 


Line aR 1 pen 


B i-es 4-37 +5} et 
and is finite for all given positive values of z, however small 
В may be, tending to the finite limit 2—1 when В is inde- 
finitely diminished. 
Let K be its greatest numerical value between 
В=0 and В=е. 
Then the portion of the integral 1 between 0 and e does not 


exceed К Ion dB, 4e. К cs and therefore vanishes in 
0 
the limit when т is at Png: increased. 
B 7 
Hence | e^ É B^ Е 1523) ap vanishes when n is made 


infinite, for all positive finite values of z. 
927. То convert 


E log '(z) e( — 1) [ 1+ tg Grp gm ‚аа |, 


the right-hand side may be written by theorem (а), 


=(— Dy» | [egre Br e- C28 pr34- M. | ав ; 
2 Фог) ав (n42) (B) 
and this includes the case 


ie logT («)= IK REM af bar (C) 


928. The same method of treatment will apply in many 
other cases. 
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Thus the sum 
Па. (p-1 
"тЫ. В?-1(е-»--е-?в--е-%е--...]@8 
E Um tig 
0 


j —¢-? 


æ B © В 
RO E OU ri UIE Ge! a.) 
nal NA x sinh 2 


i L1. 
=p gt a mT оет) 


| P" l(e- 8--e- 384 g-58.L. .. .)dB 


OS ig, ea 1 B»-1 
ape в dB = zr], sinh BOO сз (E) 


Similarly 


ЕЯ ВЕ И | 
8p—1»— apt 5p wt 


1 8? 1672 есу 1 В? -1 
=r), l+e-¥ ав= т, е? че Ыы reddi. iy 


And whenever such series occur the conversion to a definite 
integral form follows at once. For instance, in the expansion 


(Diff. Cale., Art. 574) 


sec х+{ tan х=1-- A A RW ж 4 


11! 33| 


КОЛЛ aj) | 


,-2y "P ete- Be- 38--e- 58-Le-78 | ., .]48, " odd, 
0 
Ар =з() um 8—e-798--e-9P—e-75--...]dB, n even; 
3 0 
4 nti? Bre- B 
ESI = ч pAs IC АРИ, gal н РИС (G) 
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Thus the »*^ Bernoullian number 


aon 9% 2n—1 
By = = онон р) 42-1 or yen, sili B dB; ...(H) 


and the n™ Eulerian number 


Вы (| w SEE pd. eee (I) 


T cosh 


If we write Banı аз 
айг. 1,1 ia E 
В, (2т)?" Е Но gm „. |=2@2n)! > Orry" 
=4n{ 8?n-1 (e -?*8 4 g— 4h + e-6*8-- asd) ав, 
0 


we have 


В» 4n| É зы дын" СА TO RP ccu (J) 


1-е о sinh 2898 


а, result due to Plana. (Мет. de UE de Turin, 1820.)* 


930. mt Method of obtaining Expressions for log I (æ), 


iog Г(ж), 2 yog T (æ), ... =—log T (x) as Definite Integrals is 


as follows: 


d^ 
' da^ 


Differentiating the equation г | e€-7*a*-! da, we have 
0 


me. Фое Мода decise. lud (1) 
0 
But | ede | -= | sE 
0 a 0 a 
and integrating this between limits 1 and a with regard to a, 
© еее 
loga- | E de sss (2) 


. dT (z) =| ep —Ó dà 
' ағ 0 в 2 


(0 po 
е-“—:— e—a(!+2) 
-| | etie er ERES 


0 


* See Boole, Fin. Diff., p. 110. 
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and changing the order of integration, 


pe ae Se pas i By lt al 
| 5 ө : deda=T\(2)) He рар) 


^ ‚ dlogT'(z) _ ауе 
dx 


2, 1 
D(a) дк “He 5р) © m (3) 


Integrating this with regard to x between limits z— 1 and 
=, 


D od ic dtt Da: 
ке rej- | [еше EHE | de. ....(4) 
Putting 2—2, 
ut f eu ЛАТА! 
=f е EL d 


Multiply this by z—1 and subtract from equation (4); 
© 11. 
gri) | {e-na 000 eS 
0 


2 log(14-2) 
(5) 
Now put 1-2=е’, 
mn -.g € елү dB 
log г). Салм a р... (6) 
Differentiating this with regard to z, 
d а. bul $8 
Flog Г(ж)= f. eo) DH dt (7) 
and a further differentiation with regard to x gives 
4? 
72105 Г(а)= р us „ав. . Dui aii 
Differentiating (8) n—2 times with ib to x, we get 
Z logT(z)— cf 5 ао (9) 


Results (6), (7), (8), (9) give ы, T(z), and its differential 
coefficients expressed as definite integrals. 
From (9), expanding (1—е-#)-1, we have 


T log T (2) = te wr BM e 98-- e 6*D2-Le- C22. ...) d 


1 


=(—1)” Г(”) | = Те" D tpt . to ©], 


the formula of Art. 911 lica 


www.rcin.org. pl 


| 116 CHAPTER XXIV. 


And so far as formulae (7), (S) and (9) are concerned, these 
definite integral forms are the same as those obtained in 
Arts. 925 to 927 from the result of Art. 911 (6). 


931. Approximate Summation. Maclaurin’s Formula. 
As we are dealing with many series of the form 


bius ot 
pratt: (p> 1), 


and other forms in which in some cases an exact summation 

has not been effected, it is desirable to explain the method 

usually adopted for approximate evaluation of such summations. 
Defining the symbols 2, A as in Differential Calculus, 

Art. 550, viz. such that 

Ew,—u,,, and Au,=u,,,—u,=Hu,—u, or (E—l)w,; 

and also remembering the symbolical form of Taylor's theorem, 


d 
eu, =Uz+n, Where D=; ; 


we have the following identity of operators: 

Eze], 
and it was pointed out in the Differential Calculus that these 
operative symbols obey the same elementary rules of algebra 
as quantities, viz. the three fundamental rules: 


(a) the associative law, 
(b) the commutative law, 
* (c) the index law for positive integral exponents, 


with the exception that they are not commutative with regard 
to variables. Hence, bearing this exception in mind, there is 
an algebra of operators bearing formal analogy with the 
ordinary algebra of quantities, and such theorems as the 
binomial, multinomial or exponential expansions hold. 

Let us define another symbol, У, to be such that 


LU = ии ad ug... Ua, 
where ù, is some’fixed term of the series. 
Then Уи Ўт = и, 
4.e. Ў AU, = Ug, 


and therefore Х represents the inverse of the operation A, 


www.rcin.org.pl 


MACLAURIN’S SUMMATION FORMULA, 117 


which may be written as or A7; and since А { (2) +С), where 


C is à constant and f(x) is any function of z, is equal to 
[f(@+1)+ C] — [f (x) + C] —f (4-1) — (2), 
so that the constant disappears, so in reversing the process, if 
such reversal be possible, we must restore the constant, so that 
we shall regard Zw, as A-w,4-C where C is an arbitrary 
constant to be determined in each special case. 
In this respect the symbol of finite summation, or integration, 


X behaves exactly as the sign Jas of the integral calculus. 


ш 1 pit 1 

Thus ХС Mf hrs Ma. 

Now it has been shown that 

t ё. B B B ; 

Е ta te (Dif. Calce., Art. 148); 
whence dividing out by ¢ and writing D in place of $ we have 
the following equivalence of operators, viz. 

h DIOUEEUMEDS E 

s-17p ata 274 te 

in which all the operations on the right side represent direct 
differentiations except the first, which represents an integration. 

Applying this to any function of z, viz. uz, 

he 1 В, du, В, Фи. , B,d*u, 
For this and many other formulae derived from the same 


principles, the student may consult Boole, Finite Differences, 
p. 89, ete. 


932. Apply this theorem to the case of the series 


D5—..., 


Lie 
Here — w-l, Zu-lHbtpbekÓ 
Hence 
LH ems tO [Mta as) an" 
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The constant C must be determined in such examples, either 
by reference to some known case of the summation, or by 
absolute calculation of the result for a particular value of z, 
&nd when once found, the formula can be used with the deter- 
mined constant for summation for other values of z. 

In the present case, putting s=% , 


qu. um 1 , EN 
Cp T (+++ TM +7 —log 2.) = Euler s constant — y. 


If this be available (see Art. 897) the series can be used for 
the calculation of the harmonic series to any degree of approxi- 
mation required. If C be not available take the case х=10, 
and insert the values of Bernoulli's coefficients, viz. 

B,=}, В,=з3, В,=0, B,—34y, Ву=еь, etc. (see Art. 879). 

Now 

ТФА = 2928 968 254... 

Also log, 10 —2:302 585 09; 

-. 2:928 968 25... —2:302 585 09... 


-Ü goi Gt 987389 3 $9 d$ 7 
"626 383 16... — C-1-:049 167 496; 
2. С='577Т 215 66... (Euler's constant), 
which is correct to eight places of decimals. 
Hence to the same degree of approximation we may now 
proceed to sum the series to any other number of terms by the 
result 


dij 
az 
It will be noted that to obtain eight decimal places of Euler's 


constant only three of the terms on the right-hand side affected 
the result. 


933. Take the case 
1 1 1 
atantz 


LIT jal -L— 57721506... --log,24- mi pP ate Ww 


1 1 
tint Кя (n> 1). 


1 


Неге и» 


WWW.rcin.org.pl 


APPROXIMATE SUMMATION. 119 


MAS dr 11 В, ал Be (1) ва) 
Хин я+0+ | 571-5 gi гере ын Tia z^ Y Pm 
Пр А 21 » 1 ола li _ 
~~ gt п 291 951 19 тн” 790 P 
except in the case п=1, when log replaces — LE x 
1 
— Latte A (n >1) 
dni dsl Д Пл „1. 
ae р xitg yn 12 mat 720 fn etc., 


and this series can be calculated to any degree of approximation when С 
has been found. 

In the case when a is even, the exact sums for an infinite number of 
terms are known for the earlier values of n. "The values for 1 —2, 4, 6, 8, 10 
are given in Art. 879. 

When this is the case the exact value of C is known, eg. if n—2, 


сше 
mia (Euler), and 


dive DE ud ree gend ou. MENDES I 
a6 rtm 6:3 '3035 43 1 


If n=4, C= T (Euler) and for even values of n higher than 10, 


C can be found from с={ Bon-1. (See Art. 879.) 


934. For odd indices we proceed as in Art. 932, and the 
value of the constant is to be calculated, as it is not available 
otherwise. 


Thus, if »=3, 
ati a 13 25 te 


Take the case z—10. It will be found to give C —1:202056903... to 
the first nine places of decimals, and to that approximation with this 
value of C the formula can be used for finding the sum of any other 
number of terms. 

The value of C is the sum to infinity, in all these examples, viz. 


Ж DN except when ®=1, a case which has been considered. 
r=1 
935. Consider finally the case 

log 1 +log 2+log 3+... +105 x. 


Here u,=log x ; 
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E log (z 1) - C-- log e+ | log zdz- 5log v+ 
Aug ЖЕТ 

ЖЧ cet Р. 

=O +log 2 +2 (008 2-1) ор +75 ML TN - 


1 ыл. 2 Ja d 
=C-2x+x log x+ log x+75 x 360 z+ 1960 x E 


and when z'is made very large 
log (Sez ate) =O+elog 2+ log; 
^ C=logn2r ; 


PEL jJ d I isi 
. log(1.2.3.. _ 2) = log 2m - 2+ (2+1 Jlogz-- 25 ВВ 2 toS XH 
1 


1 1 


$.е. 1.2 = ma ате е е ии QUE а, 


; 1 139 ] 
her.” 1,8 78... T Nine 1 а 5184053 ``“ 


аз а close approximation. (Cf. Arts. 877, 884.) 


936. It will be seen that the formula 
1 By du, 
и, C- [usd 0—94 Ч 
will be of the greatest service when Bd of exact summa- 
tion fail. The student should, however, test the formula for 
himself in cases with known results, such as 


— etc. 


2 2 
mud 
5 gain familiarity with it. 

Enough has been said to show that the summations we 
require in the present chapter, such as 


Mri qo NEC 1 
ани (r1 


can be readily calculated, when wanted, to any degree of 
approximation which may be required, without the labour of 
caleulating out each term separately, except for a few terms to 
determine the value of the constant. We have, for finding C, 
chosen 10 terms for the obvious reason that the arithmetical 
calculations of the right-hand member of the equality are 
thereby much simplified. 


* See De Morgan, Differential Calculus, p. 312. 
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937. A Theorem due to Cauchy. 
It is a well-known theorem in trigonometry that 


TES, 


where E, is & quantity which may be made as small as we 
please by taking т large enough (see Hobson, T'rigonometry, 
Art. 293). This is so whether 2 is real or complex. Also, 
when o» is indefinitely increased the series is absolutely con- 
vergent for all values of 2, with the, exception of such as are 
expressed by z= ул for integral values of r. 


= 22 
: т?т?— 2? 


cot Pac 
z 


Writing 3 in place of z, we have 
1 dd ад. 22 > 
ш. ЕЗД? parate 


where R'm, like Rm, can be made indefinitely small by increas- 
ing т without limit, and 


2 coth JJ ( 225) 


2 2\e—1 
and сап be written either as 
к зик. ог аз „АЙЕ. 4.e buo 
e&—1 2 ег—1 2 l—e~ 2 
1 Lo 
Hence X ris i = 
or г „аш. 
1—e7* 2 2} 
Now, by division, 
1 am AW. pus 23% 
atta at att ae gee ent 
а? 


where C EDS and is а positive proper fraction for all real 


values of 2, and the series would be convergent, and could be 
continued to infinity, provided z «a if real, or mod. z<a if 2 
be complex. 

Write in this identity a=27, 4r, бт... 2тт successively, 
and indicate by suffixes 1, 2, 3,... , the corresponding values of 
e, and let S," denote 

KC uus 1 
p'srtatedt ue 


т" 
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Then we arrive at m equations of the type 
1 1 z as g2n-2 ^ gin 
бек rey Braye CD aeg О aeg en tn 
and, adding these equations together, 
Ж 1 feo mp г ee? рае", 
Пора бт Gry КЕСУ" aee t eye n 


where ` є 342177 У) ая 73? 
and if у be the greatest of the quantities e, єз, ..., 
m 1 
, 5 | . , 
€ пни Ў нтв 4.6. є < п, 


and therefore є is also, like є, єз, єз, etc., а positive proper 
fraction. 

We thus have, taking е? to have its principal value, 
( 1 1 D- 28," 28," 28," 


7—18 2/7 бей" Ga)? * Gah 
252 


TC Ge yn gin-1.- (— D uote] R m 


and if we increase m without limit, the series S,™, S,", S,™, 
being all convergent, 


ба. to o T. T and Lt B,--0. 


Hence 
Е Wa eee. NC EHE. CEN 
(= 2-:)= Or Ga) tei : 
28, 

T(—1)-1 gem (— 1" урна". 
where O is а positive proper fraction; or, what is the same 
thing, (5—2) the same expression. 

And if we write d for ie we have 
ii 1 1 2) 
2 Ne —1 2 Z B B, B B. 
or =a p’ z?+ eee +(—1)"- "Qn! gin-2 
ey 20904 B 
А —— — А 2n4l | „әп 
i-us 2 2) +( 1) (2n 4-2)? Ө, 


where 0 «O0 «1 for all real values of 2. 
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938. Now Cauchy has shown that Maclaurin's Theorem 
for the expansion of a continuous function of a, viz. F(«), for 
the case of a real variable, still holds for a complex variable 
which is such that its modulus has a value lower than that 
for which F(x) ceases i be жее or continuous (see Art. 1299). 


The function =~ үз =—- * only becomes infinite for values 


of 2 which are m ш pius where Л is a positive or 
negative integer other than zero. "This function is therefore 
capable of expansion by Maclaurin's Theorem in a convergent 
series within the circle of convergence of radius 27 for any real 
or complex value of z, whose modulus is «27, and the form of 
that expansion has been given in Diff. Calc., Art. 148, as 

(n 3-1)- -2 5 В, ^+ Вед. .. to infinity 
ог poh mac ее, 
and the various coefficients were defined as Bernoulli's numbers. 

This series then is convergent when z is a real variable 
which lies between — 27 and +27, exclusive. It is also true 
and convergent when z is a complex variable and z lies within 
a circle of convergence of radius 27. 

And when the infinite series is not convergent, i.e. when z does 
not lie between the limits specified, the series may be stopped at 
any term (— I)7-! eui 2?^-? and the error is then numerically 
less than the next term, (— 1) "eris gin, 

This theorem is due to Cauchy. 


939. Lemma. As a preliminary to what follows me may 


remark that such an integral as [s е — dæ, where 0<0<1, lies 


intermediate between UE =p ue and IE 2202, where 0, and 6, 


are the greatest and least values of Ө between «=a and а=. 
Therefore "P йл=Ө | for some value of O between 0, 


and 6,, and therefore, if 0, and 0, are positive proper fractions, 
so also must O be a positive proper fraction. 
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940. Now we have established the equation 


ү (a) =F log D(a) = ү on “dB (Ать.980, 8); 


от, what is the same thing, 


Ve De glor cet) [ 5 ag- [ en Pag. 


‚ the finite series established by 


Hence, vali; ete for Re 


Cauchy ie 937), 


W@+1)= log r@+1)= =], 1-62 = Ps gue. 
Tope QE omiiia 
(0 — O « 1), 


al (2n4- 1) 
. — —— — — = n—127- 2n-] Bon 1 
9 E t3 eS cd @° Tto (2n)! au 


+ (— 1) By I\(2n+3) 
(2n--2)! ars 


O< Ө ЕТ), 
i.e. 


Vet I= ов (2-1) 1— UE SE. - m 


Жш, jo 
EO xe qm ms ©, (0<0<1). 


Integrating this sicht, 
me г т" 1570, 
Basis 
(2n 4- 2)2?"*? 


where 00, «1, by the lemma of the last article, А being а 
constant to be determined. 
Let z become infinite. Then 


i Г(24-1) № I"(z4-1) 
A — Lt; а) 8? |= ene та) 022] 


—(— 1)"-1 В aon (— 1)" ©, 


TR log (1+7)=0, by Art. 911 (3). 
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Hence 
d Ic Ф 
(e+ 1) 5. log Г(®-Е1)=1од 2+5. CIN ^ш 
n— Bony n Bon +1 
—(—1) 'Sagin (ПО n 2)aio Ov 
(0 « 0, « 1). 
Again integrating, 
, Bol 
log l'(z4- 1) — A'4-z (log x—1)-- } log 2-25 т 2 
Boa By 1 


+(-1)! ө, 


Qn ^j aie C Den p zen 


(0<Ө,< 1), by the lemma, where A’ is a constant to be 
determined. 
Let z become an infinite integer, 


A' — Lt... [log I'(z-- 1) —z(log z—1)— 3 log x] 
= Lt... [log (/ 2x7 x*e-7) — (x -- 3) log x4- 2] 


—log /27. 
Hence 
log l'(z4-1)—1 log 27 4- (z-- 1) log z —«-- —s B, 1 rif eire 
1.22 83.423 
a1, Вы n Boa 1 
TOOL н) ан акт (2) би) (2n-E2) ӨР 
(0<Ө,< 1). 


This result is also due to Cauchy. 


941. The series, if carried to infinity, is known as Stirling's 
Series. It is divergent, however great z may be. For the 
general term 

Bai 1 1 1 2(2n)! 
(2n—1)2n af (2n —1)2n gai (дт)?" 2m 
and the ratio of this term to the preceding term is 
(2n—3)(2n—9) | Son 
(272)? ge зи 


2 
i.e. ultimately E and however great z may be, will 


ultimately be > 1 when % is large enough. "The formula can, 
nevertheless, be made useful for approximative purposes for 
calculating I'(z4-1) For, as in the series of Art. 938, the 
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error in stopping at the term involving has been shown 


gina 
Baa 

to be Ө jrm (0< ә « 1) i.e. the error is less 

than the succeeding term. And as the ratio of two con- 


(2n—3)(2n—2) San 


secutive terms, viz. Ora) быр is less than unity 


until (2n— 3)(2n—2) s Son exceeds 47?2?, the absolute values 
2n-2 


of the several terms go on diminishing until this happens, and 
then increase again. Hence the closest approximation will 
be obtained by continuing the series until that term is reached 
which precedes the smallest term. 


942. We have as successive approximations 
log l'(z-- 1) > 1 log 27 (2+ 3) log x—z, 


log I'(z4-1) < 3 log 24+ (4-1) log a— as Б 
(SUP WE 
log '(z4-1) >} log 27 -- (z--3) logz—2--175 73-4 zw 
log l'(z4-1) < 1log2z--(x--1)logz—2 
21 Be I Ld 
*pàz 8.447 5.62» 9 
1 
And since B, =в› p 2,30. B, gg ete. 
IXz4-1) 
> J/Qrxx%e-®, 
1 
< У/2тхх®е-=еї, 
1 1 
> J/2Qraarte-tele 360, ete., 
i.e. 
l'(z4-1) 
>N rar e”, 
LL ара M i 1 
«insta (газаа), 
ыч Lr ME: 
> ттан (1+ tp mu i oben aad 


ete. 
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943. In order to facilitate calculation from the series 


log I'(z4- 1)- log 2m 2+5) log 2—2 


it is desirable to arrange so that ж shall not be small. 
For this purpose Legendre puts х=4--а; whence 
log '(z4- 1) —log x+log I'(x)=log x 
+log I'(a)--log a(a+1)(4+2)(a+3) 


and 
1 1 Bi dil ЖР 
BNI i 
ES у we ЮБьа(а-+ 1)(a+2)(a+8), 


where д is the modulus of the logarithm tables, viz. 
w= logype= "4342944819... . 

Thus, if log,, T (1:25) be required, 2= 5:25, and 

log T (125)—5 log 2-475 logi 525— 4525-25 z5 
—logo[(:25)(2:25)(3:25)(425)], 

and,by this artifice it is possible to avoid the calculation of 
all but the earlier terms of the series We could make 
х=б-ра, 6+a,..., equally well, and the choice is in the 
hands of the calculator. 

Legendre remarks as to his calculations of the seven-figure 
tables of log Г (х) with regard to the above: “de cette maniere 
on n'a jamais eu besoin de caleuler plus de deux ou trois termes 


2 Ne M. 6 " 
de la série 1.29 4$ 3.4 5 ee et pour avoir logT (a) 


— etc. 


approché jusqu'à sept décimales, dans tout l'intervalle depuis 
а=1 jusqu'à a=2” (Exercices, p. 300). 

Legendre's m, К, А’, В’, C' are what we have called џи, x 
B,, В,, В, respectively. 


› 


94+. The Case when x is а Commensurable Number. 
We have established the result 
p 


© log г) |" (5 - лай dB. (Art. 930 (7)) 


=e 
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And we have seen that Euler's constant y is the value of 
— Flog T(x) when х=1 (Art. 911 (4)) 
А (=: P drach 
that is y--[ C5 mee 73) dp. 
Hence, adding 
d е-В—е- вх 
qe OS Г(®)+у= : ree 4B. 
In the ease when z is a commensurable number* this 
integral can be reduced to the integration of a rational 
integral algebraie expression, and the integration effected in 


finite terms in terms of the ordinary algebraic, logarithmic 
and inverse circular functions. 


Let = ‚ Where p and д are positive integers, and let e-8— t4. 


loda. fp 


d 
Then z log I'(x)+ у=, Tim 
and the integrand is a rational integral algebraic function of t. 
If q=1, i.e. if x be an integer, the value of Ž log Г(2) is 
given by 


d 11—12-1 
215198 re+y=f qe. dt 


=f (14-4 ... +t) dt 
0 


EDIT 
илаа duc" 
as might be expected from Art. 911 (2). 


945. Expansion of Г(ж-Е1) derived from the Integral Definition 
(De Morgan). 
The expansion of log Г(1--2) in powers of z may be obtained 


directly from the definition of l'(14-2) as | е) dv. 
0 


For we have ты . T) eut. 
— 9-202 
Hence T'(14-2)— Lt, ШЕ PIN E. 


* See Serret, Calc. Intégral, p. 184. 
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Let e-^—y. Then adv= -%, and 


Г(1+2)= Lt al, ge eb 


= [4 PN. [> 2+ 1) (Let Ih, a positive integer.) 
— Lis be ГОГ), 
Г(0+2+-1) ' 
r(b. 
l'(z4-5- I) 
‚и @+®@+Ь—1)... (E D T +) _ 
a nd (b—1)(b—2)... 1. Ь=+1 


i.e. Lt. boc =1; 


1.6. 

2 z : 
log '(14-z) [4 E log b—log (14%) —1юв (1+2) Бы cd 
or, expanding the Нч assuming 2 « 1, 


log (14-2) =| —({ i +5 ер +5 logb)z 


tanta +p eg (ptt tg) t] 
and when b is indefinitely increased 
2% x g^ 
log Г(1+=)= —}у®+5»%- — 8, g tS, qer 


for values of z, 0<e<l. 

This investigation is due to De Morgan.* 

It was felt desirable to deduce this series directly from the 
integral, rather than to base it upon results deduced from the 
property I (x+ 1)=xT (x), i.e. the difference equation t, —2zu,, 
inasmuch as Legendre’s tables of the values of the Gamma 
function are derived from this series and others obtained from 
it. And in default of direct derivation of the series from 
the integral itself, some doubt might be felt as to whether 
Legendre’s tabulated results were the values of the integral 
itself or the values of the integral multiplied by some periodic 
function of z whose period is unity, which, as explained in 
Art. 863, would equally be a solution of the difference equation. 


* De Morgan, Diff. Calc., p. 584. 
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946. From De Morgan's investigation given above, the 
formal identification of Г(24-1) with П(2) for all positive 
values of z, may proceed as follows: 


1 (2)=Ltyamut (155) (1 +3) ти (145; 


4 log пам, -„[з log p—log(1 +7) —log (1 +)- 


—log (1+=) |, 
S, 


and if z« 1, aa К A j 2% Sa... 


'. II(z) 2 T (z4- 1) if z «1 and positive. 

If z lies between 1 and 2, say z—14-£, then, since 
Ws ыы and II(£) -T(1--£) (0 £« 1), 
it follows that II(12- £) 2T (24-£), 

1.е. П (2) - (14-2) when z lies between 1 and 2. 
Similarly if z lies between 2 and 3, etc. 


Hence, for all positive values of x, II(z) and Г(1-+2) are 
identical. 


a 
947. The Integration zl e€-*v" dv, (a not infinite, n> — 1). 
0 
In considering the-integration of e-'v^dv between limits 
О and a, where o is not infinite, we must have recourse to either 
(1) an expression in series 


or (2) a continued fraction. 
"+1 1 а 
==! e-"v^dv-2l e-*—— "m —vant+l 
(1) ial e "фу = L ЖТ, еке ?yntldy 
е-аа%+1 1 
= ple 
and by the continued use of this rule, 
e чат аз 


а? 
"AGE ты *ix:3)(n4-3)(n1-4) 
+... аа inf. |, 


I,— 


www.rcin.org.pl 


INTEGRATION BY CONTINUED FRACTIONS. 131 


а, series which is always convergent for any finite value of a, 
but only slowly so if a be >1. A little consideration will 
show that the integral remainder is ultimately infinitely small. 
Or we may proceed thus: 


Let J,= [чао есе rem XM 
—e anc ти, ; 
whence 
PP tu tmm D, qnin De (nr) 
J4,-—e-7*a I = | 


+n (n—1)...(n—7)4, .,.,. 

If n be a positive integer, the integration can be effected т 
finite terms. But if п be negative or fractional, the series on 
the right-hand side is divergent if continued to infinity what- 
ever а may be. The terms however ultimately take alternate 
signs, and when such is the case, and when there is convergence 
for & certain number of terms, and then ultimate divergence, 
we can apply the principle adopted in Arts. 938, 941, the 
convergent part making a continual approximation to the 
arithmetical value of the function under consideration, and 
the error being less than the first term omitted.* 

If then J, be thus approximated to, 


mete f -[) e-* y^ dv, 


and I,—DI(n4-1)—7,. 


948, (2) De Morgan has shown how such an integral as | е v" dv can 
be converted into a continued fraction. у 


When this is done Г е" = Г(п+1) -Í e ^" v" dv, as before. 
о v 
Let [ew "dv—e-*v" V, where V is some function of v. 
Then differentiating with regard to v, 
—e-*y"=e-"y" Р'+ ne~" y"! Dies e-*y" V ; 
~ vYV'- nV -vV = —v, 

or vV'—(v—n)V — v. 

Consider the equation 

OVS ЕЁ, uult rte (1) 


*De Morgan, Differential Calculus, p. 226 and p. 590. 
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Putting V= , we derive an equation 


1 +h, Жи 
vV; —(v—a;) У, vtd 1,1, уе if оа: e (2) 
where b, —a,=k, b, =k, =, – а, а= — (a4 +1). 
Putting уі in equation (2), we derive an equation 
1+6? 
vV; —(v—a;) V,—v+b,V,2, 0060000606070 66000000006 (3) 
where b,—-a,—k, b,—k, аз= – (а,+1), 
and so on. 
‚_ 1 kv ky ky 
Then lero te th Te Oa. 
In our case 
а =", biz:0, k= -п=0, ; 
a= —(14-n), b= –т, Ё,=1=6,; 
az =N, b,—1, k;= —(n—-1)=& ; 
= —(n+1), b= —(n- 1), ky=2=bs ; 
a,=n, b,=2, k,=2-—n=),; 
etc. ; 
whence 
pablo “Аг 1 nv v^? (n—1)v 2071 (n- 2)v7 (n— е 
ета [т А 


The expression converges rapidly for large values of v. 
The process above employed by De Morgan is similar to that employed 
by Boole, Differential Equations, p. 92, in the solution of Riccati's equation 
d 
2 Z — ay +b’ = 
The equation we have just*solved is a very similar equation, viz. 


d 
22 kay - by? = – 2+0. 


949. More generally, consider the differential equation 


P+Qy+ Ву+89 =0, 


ах 
where P, Q, R, S are functions of x alone. 
Let X,=Azx*, Х,=В28, X,=Ox", etc. 
Take y1, Уз, Уз, ... Successive new dependent variables, such that 


Xi Xs X; 
Ir! УТ у Ту, 
Then when A, B, С,...а, В, у, ... have been properly determined, we have 

Аз“ Ваё Оту 
71x TF TH.. 
viz. a solution in the form of a continued й. [Lacrorx, t. II., р. 288.] 


у= ‚ etc. 
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To begin with, using accents for differentiations, 
X61 35)- Xi, 
(AX i 
. PIL X; +В Хү X'403)-Xii. 


ї+и ^ (ry? QA» 
ie. (Р+090Х,+ВХг+8Х’,)+(2Р+9Х,+8Х',) у. + Py? - ВХ у, =0, 
or Put Quit Ry +81y'1=9, 
where P, P+QX,4+ RX+8X’,, 

Q,=2P+QX, +8X’,,. 

R,= P, 

8, = — SX, 


At the second substitution, viz. y,— 
becomes 


i the differential equation 
Pj Qin Roy? +S2y'2=0, 

where Р,, Qz, Rz, S, are formed from P,, Qı, R,, S, in the same way as 

the latter were formed from P, Q, R, S, and so on. 

Again assuming the expansion of y in powers of «x to be of the form 
Аха -- A,ze*l-F... and the expansion of y; to be Bz8-- Byx8t+1+..., and 
so on, we can by substitution in the several differential equations they 
satisfy obtain the values of А and a, B and f, etc., by an examination of 
the lowest order terms occurring, and thus express y in the form of a 
continued fraction. 


950. Development of Y (a+ 2») =“ log (a+) in a Factorial 
Series. 
Since 
Ау (a2) m Vr (a- 4-1) - (a2) =: flog T (a 2-1) - log Г(а+#)] 


d 1 
= 48.05 (@+4) =, 
we have 


1 1 (= 
A?y (a +2)= А а+ж+1 atx (а+х)(а+х+1)”; 


i PEE P TTET DN 
Ачаа) s аар ET УЗу 


and generally 


A y (a 4- x) = Д" 


Let 
ж? gm 
y (at x)— A ЖУ oat Al is Ron ЖА 7+. , 


(- 1t (8-1)! 
rz ~ (ac z)(ad-c-1)...(a-z4n-1)' 


where 5” zz(r—1)...(x—-4 1). 
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gin m 
Then Ay (az) — A, 4, 11+4:51+ 
- S 
Д?у(а+=)= А+, П TP + 1 
etc. 
Hence 


A,— (a0) A =Ау(а+0),  A,— A*j (aO), ... etc., 
where A"y(a +0) means the value of Ду (а +2) when z is put =0. 


Hence 


аа = d log Гана - (a) ет 


 lz(z-1)(z-2)(z-3) 

4 a(a--1)(a 4 2)(a 4-3) 

a series which will terminate in the case when z is a positive integer 
and is in any case convergent for real and positive values of z and a. 


+..., 


The value of у(а), i.e. S log.T'(a), can be found for any particular 


value of a by means of the series 


d ly 8,1. B 
ы, == es eis ж тш. дыг 
Ta 108 (c 1) 105,2 2+5. T 1 etc. 
of Art. 940. 
951. In the case when a=1, we have 
H i.d —1) (2—2 
уча) 1 2—1) 1 е9) 
_1=(#-1(=-3) (2-3), 
4 4! 
and — V (D) = y (Euler's constant). 


Since Az = nz'^-", this may be written symbolically as 
1 1 1 1 
и yt atga- sa) = - y * A log (1 +k) %, 
к а Е 
i.e. 25198 Г(1+2) = -у+А log (^)= 


952. Other properties of the y function are: 
Since I' (x - 1) = Г (х), we have by logarithmic differentiation 


WEAN) = mi. ette (a) 
Since r(x) r (1 -2)=—"_, we have similarly 
y (x) — V (1—2)2 #4 COL zs. ................. ee (b) 
Since 2?* T(z) T ($+2)=2\r I (2x), we have similarly 
V(r)-V(à42)-2V(22) -2log 2. ..................... (c) 
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Jr (x M 
3 )-Zr H we have similarly 
cos = 


v(x) 4 (155) - (2) 10g 245 tame ............... (d) 


Since T (2) T (+ +т)г (2+2) Lip (a 1) аннат) (na) 
we have similarly 


уа (2+1 +++ (++2 +. +0 (24+ —*)= ny (nz) -nlogm. (e) 


Все ача r( : 


953. The equation дуаа) = 1 is of considerable service in 
summation of series. 


l. A sum of the form 


1 1 


1 i 
гы at aa to n terms, viz. 


S= X —— can be written 
1 


2. A sum of the form 
] 1 1 1 

one. ak 3b * a 3b 3b are M 
UY he У 


ste) арын) 
=g" x 218 = 


Eg. о 


УЗИ: Eavaio-i[vasn iva m-vox 
0) 3- H3 cad inf 
= ip i ire 
ауду) дув +) 
=i[va+n] -#[%@+5] 
Муф -уФ) 
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But by (0) (==) У(Ф-У@=т; 
-. the series is-7 А 
which is well known otherwise, being Gregory's series for tan^! 1. 


3. Sum the series 


Mi db o qa 
Here ; 8-19 .- iX У D 
-)EAV +7) FEA +r) 
=i[va+n] -3[va so] 


-41UQ0)- 0). 
Now by (с) (#=3), — V(0)*V0)-2/(1)-21og2 ; 
n V(1)-V)-21og2 ; 
. S=log 2, which is well known otherwise. 
We may note that it follows that 
V) = V(1)- 21og 2— – y – 21052 
= —0:5772157 — 1:3862944 
= —1°9635101.... 
Ву (с), УФ + ()=2 (3) -21og 2=2{Y(1) - 2 log 2) – 2 log 2 
= — 2y – 61052 
and VQ) - 0) =т. 


Непсе ($) = $-Y-3log 2, 


уф = -$-у-31од2 
and y (4)= - y -21og 2. 


954. Gauss has established a remarkable result, giving for the function 
y(x) the value of (1 — 2) + y(x) in a series of trigonometric terms in the 
case when v is any commensurable proper fraction. This result taken with 

(1-4) - V (x) 97 cot er 
will enable us to calculate the value of у(х) in all such cases. 
The theorem is given by Bertrand in Art. 307 of his Calcul Intégral. 
For shortness we shall denote 


log x by £z, v() by Yp, cosr@ by cp, 105 4віп? T by Z,. 


Then when 0—2" or 4r or 6r ГӨР Bahtir 
[ q q q q 
(470. —034,—...—1; Сар 04a pm. me; dO. tHE $e, = 
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Writing the fundamental équation 
1 
у(х)= Lu-.|logn-i-—- ра Jas 
ar n 1 n+l 1 
ve)=(-2) 4 (19-255) QA) n Зы 


and putting == where r № 4, and both are positive integers, we have 


„(а -2)+ (e)k. + ыы S 


Кут r=1, 2, 3,...9 in this equation, multiplying by cos б, cos 20, 
cos 30, ... cos 40 respectively, and EOR we get 
g 1 
NOM -H с, «($ 05) рч 2с). (Хе, a Exe) а 
Now the coefficients of log 1, log $, log 2, etc., all vanish and since c,=c,,,, 
etc., the remaining terms form a continuous series to infinity, viz. 


ELS UE 


^ Хер, =i А , 
1 
viz. an equation connecting y, V», Vs, ... Ya-1s We, the last of which terms 


is у (2) о) — y, where у is Euler's constant. That is 


CW + €3 Vra + 65 Vra +... + Cg у = Ly. 


4т 


So far 0 has stood for any of the quantities zi qt or 92-1 7 т. Say 


the first. Then similar results will hold for the rest, ùe. if we take 20, 
30,...(q—1)0 in place of 6. We thus get q—1 linear equations from 


| г) (3 пит) 
which we can find (2) (2) Ses »( 77 viz. 


Crit са +...+ CpPpt...F Gp Var +... + ааа 1+7, 


A 


Cait саа+...+ Capp Cog Yop +... + ско 7 Ea y, 


Oa Cea. Capp. Ca p Pap tee салаа slaty, 


озона ооо вонь ооо оное но errr err а ооо ООО но оао ое о ооо ооо о оо свое о оо оо er ОТОН 


64V + ‚++ беләр}... C e plU ap + + + Cg 19V o1 = gai ty; 
and in addition we have 
(UA T Coa... + орф... Наиль + +++ + Саад) Мал = —(4- 1)y - qlogq, 
which is merely a case of the identity (e) of Art. 952, for the coefficients 
cos 90, cos 240, etc., each =1. 
To solve these equations we multiply them, and the identity, respec- 
tively by ср, Cap) Caps +++ Cops 
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Now note that cc, + сол €), C3, Cg, +... + со Са for any integral values 
of А, р (the last term being unity, since 40 =а multiple of 27) 


=} Eesti È caur ; 


and that each of these sums is zero, except in the two cases А+д=а 

multiple of g, and that in the cases we have to consider А and p each range 

in value from 0 to 9-1. Hence the only cases of this kind are when А=д 

or À—q-— р, and both would happen if A—4-—q- p, i.e. if g be even, and 
MAR LT. 

A=p= : 

If A=n, Ўсе 51-1 ; ifA-g-p, tÈ casur =+21=$, 


and when 4 is even and л=и=1, Eco Е иу. 


The latter case will occur when, q being even and therefore q— 1 odd, 
there is a middle term in the system of unknowns, viz. y, —, ,— (3), 
and the case need not be distinguished from the others. ‘Thus, after multi- 
plication by ср, сз, ... Ср and addition, the coefficients of all the unknowns 
vanish except those of V, and p,p, and the coefficients of these terms аге 
each 1 ; and if q— 1 be odd and p=, all vanish except that of (3), which 
is the middle unknown of the series, and the coefficient of this term will 
be q. 

And on the right-hand side we have 


(o Lit eL... + C4 p Гал) + y(Cp + Cap + .-. + Cap) — Ven – d 084. Сар 


=F (cp Ly + cap La +... + саль Гаа) — Vy - 41084. 


In the bracket, terms equidistant from the ends pair, but if 4 be even 
there will be an unpaired term left in the middle of the series. This term 


is { соз 2р6 log 4 sin? g which reduces, since q8— 2r, to q( — 1} log 2. 


Hence tbe right-hand side becomes 


е 


(0 +сь +... +з, ал) ау 41084 (9 оаа), 


ог a(o Ly сз» 11+... +з, а-а) ay —91089+9(— 1)? 1059 (q even). 
2 2 
We thus have 
(o2 
y(1 t. v(2)=2 { > ep L,- Y- logg} (4 odd), 


(32 
or =2 у ees Lo — y - log (7 1)*log 2} (4 even), 
1 


and this, as pointed out above with 


v( -2)-y(2)=7 cot? т, 
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will enable us by addition and subtraction to obtain both 


-2) and (2) 
SILEAS 
for any integral values of p and g (p< q). 


It will be observed that these theorems give the tangents of the 


slopes of the curve y=log T (x) at equal distances on opposite sides of 
the ordinate at z—0'5. 


Ex. If p=1, q—3, 
т т 
¥(§)-¥(})=7 inui тс 


V) (3) — | -y -1og 3- cos бу log 4 sin? | 
=2[– у – 1053 – $1053] 
= – у – 31053; 
> = -у- #108327, | 
«АЙ Рр Df Me 
V= -у-#1о3- tr. 


955. List or RESULTS. 


As the results obtained in the present chapter are very 
numerous and necessarily scattered over many pages in the 
gradual development of the theory of Eulerian integrals, it 
may be convenient to the reader to have the principal facts 
arrived at collected together for ready reference. A synopsis is 
therefore added in two groups, the second group referring 


more particularly to the y function, which entails some 
repetition. 


GnouP I. 
Вет), p-| gi-X(1— zn dz. (Art. 857.) 
2. If l, m be positive integers, В (1, ран: 
If Lonly be a positive integer, 
E. (1—1)! " 
Bil, Wy ө ү рү эте лу үү чене) 


gmt d 2-1 
4. f (z—by-1(a—z)y"-! dz— (a — by*"-t B(L, т). 
d (Art, 859 (4).) 


3. Bil, m)=| dz. (Art. 859 (2).) 
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: re) (277) 
5. [зоб costó ао 2 > (Arts. 859, 869.) 
а 2r ст 1) 
^ 6102-10 cos*m-19 а 
я. го) [ arte dz, Г "ni a"—le—kx dy, 
Ж: ^ kn 0 
1 
1.2. 
vica +) ‚ П(2) = Пи 
d Im d eG FIFY -AF 
Р (Arts. 854, 864, 874, 889.) 
8. Г(љ+1) 2A T (n)—II(n). 
II(n4-1)—(n4-1)II(n). (Arts. 860, 890.) 
9. T (4)=/7=I(—}). (Arts. 864, 882.) 
10. Г(2)Г(1—2)=т cosee ar —II(—2)II (z— 1). 
_ Га) Г(1— 2) = гт cosec т. (Arts. 872, 893.) 
T gi-1 
ү? (2х) 3 P 
12. г(2)г(2)г(® Joe no) (27) * (Art. 873.) 


n 
-1 


13. af odi 2+2). (a *—)- T(nz)(27) # nè, 


Гоген ГО») reth г(2+2)- m ren. 
(Arts. 903, 905.) 


1:23:8.. 
Art. 877). 
14. Ге PIU Li 
` Pat) 207.3 
ib e He -55 3: (Art. 884.) 


16. у=0°57191566...= Lt, „(++ D +2 —log n) : 
(Arts. 897, 917.) 


ttn z T? Acl n 
17. | * о Өдөр. m | 
sis — 


(Art. 910.) 
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18. 


19. 


26. 


27. 


28. 


2 а” 
т log Г(1+2)= —yz- $5 —& 5+8, 4— 
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d mk ok 
> j; 08 Г(2)= Ltn- af log n 5 — zz | 


хн Z2 at e f. kk, 911 (8)) 


d 
=, log Г(2)= EL аа 04 inf. 
da 2 m Grm (Art. 911 (1)) 


‚ Lt, (EC) 1g n)=0. (Art. 911 (3)) 


T(n) 


(Arts. 911, 916.) 


. log '(14-2)— log ..^ =, 7 ——tanh7iz 4-(1— y)z 


sin жт 


erga. ‚ (Art. 919.) 


. Min. ordinate of y=T (x) is at х=1-4616.... (Art. 922). 


ше re-| [e-ne- 5-2 |. (Art. 930 (6).) 


1—e? | В 
i b log I'(z)-: R (2-1 on) ав; (Art. 995.) 
кув = f. fes T Е (Art. 930 (3).) 


Z log r)- c-iy[ £ «ды 48 (n42). (Art. 930 (9).) 


S75 a tmp m. d, 
ss star), sai add 
UT CR CUM iP к (Arts. 928, 929.) 
Poi (pe Dye LI 8-2], кв ® 
AAPA XM 48. (Art. 929.) 
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99, Уи =C | dz B, du, B, u B. 45 
$ 6 fu- 79 TEL | Ho tm 4| Lc Ux 


d b d SE 
30. j; o2 Alani. A 


T E 


Zin Ema То (0-021). 


ihe Sa (2n-+2) A 
(Art. 940.) 


31. log (z4-1)—4log Эт (z-- log 224-5 l P +... 


By ay < 1 
Qn—1)2n ga 1)” Gai) (n3) at? 
(00-1) (Art. 940.) 


Db. i, ТИ а TE 
J9-za*e-* | 12a (1220) 30(12z) 120(12х} 


See also No. 15. (Art. 942.) 


+1 


32. иг. 


956. II. Group ОЕ y» FORMULAE. 
Since the y-function, viz. уба) =. log Г (2), is а very 


interesting function, and very useful in itself, we gather 
together the principal results which refer to this function 
in particular. 


-re uae NEE E 145: 
1. y(z)= нь log n z--1 y" sul 
(Art. 911.) 


2. (0)——, W(l)=—y, W(14616...)=0, (ac) ә. 
f (Arts. 911 (3), 922, 923.) 


; Bep brin 
(Art. 911.) 


M 1 
4 Vm oi ips (Art. 911.) 


7 С 1 ав 
5 ўа) [e i573) ав= [| os] 
(Arts. 995, 930 (3) and (7).) 
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20 


e 


10. 


11. 


20. 
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(a= IR £t 7, dB. (Art. 930 (8).) 
: V@+1)=log 2+ Ва... (Art. 940.) 
‚ү @+1)=1- E MEI. (Art. 940.) 
\ yG--vy- |, тад (Art. 944.) 
у(2) -у= = 6 lt a (x integral) | (Art. 944.) 
d - a) E [а (From 10.) 
i Ay (a) 21... (Art. 950.) 
А уба 1) уа) =. (Art. 952.) 
. V(1—2)— y (z)—' cot zT. (Art. 952.) 
. V(3o-2)—V(43—2)-—7 tan zr. (From 14.) 
уа) уа) = (22)—2 log 2. (Art. 952.) 
уба) iv ==) = (2 )-tlog 2-45 tan, (Art. 952.) 
Wot (2) e (o3) 4-1) 
= тр (nz)—n log n. (Art. 952.) 
"xut x 12(2—1) lz(z—1)2—2) | 
j оле а а(а+1) 3 a(a+ 1)(@+2) — "m 
(Art. 950.) 
p p 
va-t) l 
9-1 
=2| y-(1)—log q+ >) cos Li al log 4 sin? =] (q odd) 
т 4 (Art. 953.) 


Т log 4 sin?” T |H- ny2og2 


2-2 
=2| 4(1)—log q- Y соз "Р 
; (q even). 
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957. Table of Values of S,= ijs sta 2 Ln. . ad inf. 


up to p=35, which is the "wis in ifta iy tenth decimal 
place is affected ; all remaining ones to this approximation may 
be regarded as —1. (De Morgan, D.C., p. 554.) | 


р S, to sixteen places of decimals. 

‚1 | 057721 56649 01532 9... +108 © (Euler's Const. + oo) 
2 | 164493 40668 48226 4 
`3 | 120205 69031 59594 3 КЕ 

4 | 108232 32337 11138 2 8T 

5 | 1:03692 77551 43370 0 ES 

6 | 1:01734 30619 84449 1 833 

7 | 100834 92773 81922 7 3 

8 | 100407 73561 97944 3 аг 

9 | 100200 83928 26082 2 Bu s: 
10 | 100099 45751 27818 0 о. 
11 | 100049 41886 04119 4 SE $8 
12 | 1:00024 60865 53308 0 =. 8% 
13 | 100012 27133 47578 5 ae gA 
14 | 1`00006 12481 35058 7 “у - T1 
15 | 1:00003 05882 36307 0 mn ФО: 
16 | 100001 59822 59408 6 BE gge 
17 | 100000 76371 97637 9 “з 888 
18 | 100000 38172 93265 0 b MP e: 
19 | 100000 19082 12716 6 Ня: Заз 
20 | 100000 09539 62033 9 $9259 
21 | 100000 04769 32986 8 SESaSS 
22 | 1:00000 02384 50502 7 a BR d 
23 | 1-00000 01192 19926 0 qagzgore 
24 | 100000 00596 08189 1 BS. Ea! 
25 | 1:00000 00298 03503 5 чая 
26 | 1`00000 00149 01554 8 oR: 86 
27 | 1`00000 00074 50711 8 2s: 2 
28 | 1:00000 00037 25334 0 вез % 
29 | 1:00000 00018 62659 7 Bs - 
30 | 100000 00009 31327 4 бое HB 
31 | 1:00000 00004 65662 9 е3" 

^| 32 | 1:00000 00002 32831 2 É add 

33 | 100000 00001 16415 5 Ze 


34 | 1:00000 00000 58207 7 
35 | 1:00000 00000 29103 8 


PROBLEMS. 
1. Show that (i) roro- di) PQ) r(5) = r (2). 
2. Show that 35 (T (2))? 2 2523 T (1). 


3. Show that P'(-1) P(2) Г(3)...Г(9) бт. 


www.rcin.org.pl 


PROBLEMS ОМ GAMMA FUNCTIONS. 145 


4. Show that 2%Г (n - 3) 21.3.5 ... (2n - 1) Az, where n is а posi- 
tive integer. [Oxrorp IJ. P., 1888.] 


5. Show that T ($ — x) Г(3+2) = (1 - 22) r sec zz, provided 
-1<22<1. 


6. Show by means of the transformation zy = 0, y = и +v, that 
TETUR um 
| ; | A arent — de dy = Bm, я). 
[Corr. y, 1901.] 
1 
7. By means of the integral | gm-1(] — z*)^ dz, prove that 
0 , 
ду: cse фиш... де А 1 бы а 
(m)n! (т+а)(%—-1)!1! (m«2a)n-92))2! ^" $ (m 4- na) т! 
а? 
© m(m +а)(т + 2a) ... (m 4- ma) 
(Sr. Joun’s, 1884.] 


!г(*) 
ar (> tnl) 


Show that this integral may be expressed as 


8. Show that the product of the series 


luin" са. di2i;3.9 1 
са Б т В + etc. 


31711514 733 9,456 149 
4 11d L5 Y bled ade 
T 9° 9 95 "4 4 ii 2,46 Bt Р 


[COLLEGES а, 1883.] 
9. Prove by the substitution 2? = £ that 
ё т 
ет = e722 72n+1 dr 
|. zonde = I LÀ T1: 5 h gent! dy, 
where n is a positive integer. 
[See also Art. 223 (5).] [CorrEGEs a, 1890.] 


10. Show that if К be any positive constant, 
K (K— 1 K 
| | ume 211—1 gm dx dy = | ( 1 РЕ e) ym dv ы | et ultm-1 du, 
0 0 


o Jo 
and by proceeding to a limit express B(/, т) in terms of Gamma 


functions. [Oxr. II. P., 1902.] 
11. Show that the sum of the series 

1 1 т(т+1) 1 т(т+1(т+2) 1 

il tae dM qud 31 nid 

is T'(n - 1) (1 - m)/T (n - m4 2), 

where n> — 1, and m < 1. [CoLr. y, 1899.] 


ө 
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1 
12. From the value in Gamma functions of | sin? 0 соз 0 d6, 
0 


show that »т(?®1)г (2+2). УтГ(р+1) 


for all real values of p. [Тьтктту, 1886.] 


13. Prove that | e—? їх =e% х 0:09811 nearly. [Trinrry, 1896.] 
5 


14. Prove that 


n 
(re) (1+3) (1+3 
лө, ... 0 © 


т n\ n n 
(1+1) (1+3) (1+3) [Oxron» II. P., 1888.] 


and Г(л+1)= П 


r-l 


1+ 


^ am" 
4 E- Кы 
>м 


[Oxron» II. P., 1903.] 


15. Show that, when z is positive, 
T(x) e 2n! 1 
Г(=2+ ll) »-02?^ n! n! z - ^ 
[Матн. TRIP., 1897.] 


222-1 B(x, 2) = Jm = 


16. Prove that, if z be positive, 
1.3 1.3.5 = T) 
м 


[ГМАтн. Твтроз, 1897.] 


17. Show that, when z is a real positive quantity not greater 
than unity, EA 1 
ь (x) = (2) + 5 „2(®+1)(®+2)... (r4 n) 


where f(x) is a function of x not greater than unity. 
(Maru. Тктроз, 1897.] 


18. Ил lie between zero and unity, prove that 


т 


а. т 
ta n = 1 — m, 
ad [CoLL. а, 1890. ] 


19. Show that the perimeter of a loop of the curve 7" = a"cos nô is 


(га) |) 


www.rcin.org.pl 


PROBLEMS ON GAMMA FUNCTIONS. 147 


20. Show that if z, y be a point on the ellipse z?/a? + 12/0“ = 1, 
and 2r be the conjugate diameter, and the integral be taken round 
the whole perimeter, then 


r4- 3 i 
IE p 2 Ж. 
Г(/+ 1) ab  [Conueces, 1892.] 


21. Express in Gamma functions 


1 Д 
| (1 —a™)» dz. 
0 


(Trinity, 1896.] 


22. Express іп Gamma functions the area of the curve ус? = а2° 
(с> 0) for positive values of z (0 to oo), also the volume generated 
by its revolution round the axis of z. [Sr. Јонм”в, 1883.] 


23. П 2sinnz Г(п) $ (n) = (27)^ $ (1 - т) ((— n 07) where 
=; – 1 and $(n) is some function of п, prove that 


n n 
г(5) 2-14 (n) 
remains unaltered when 1 — » is written for n. [CoLLEGES a, 1881.] 
24. Prove that 
< „ГАГА. 244 244% T4 
Cre ТЕ түк p ran | 9099 dua 


[Dr Morean, Diff. Cal., p. 591.] 
25. Prove that 
е oly? sci 9071 2071 Set Bet 
e~ log v dv = е" 


v 


[De MORGAN, p. ma 
26. Prove that 
lz(z- 1) ,12(z- 1)(2-2) _ 
4 L4 9. L.2.9 
[Dx Morcan, р. 593.] 


d 
2808 Г(1 +2) = -ytt- 


27. If ф(х) = is log Г(1 +2) and х be a positive integer, show that 


фб) =Ф(0) +++. +1. 


Prove further that 


$(0) = г. €-*log x da, 
and has a finite value. [I. С. S., 1898.] 
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28. И (1+2ж)”=1+А,®2+ 4,22 +..., where n is any positive 
quantity, prove that 
2% '(n 4- 3) 
2 2 mem. eeu у 
1+А4,%+ 47+... Ja Tn +1)" 
[MarH. Trrpos, 1895. ] 
29. Prove that if 


Ла) = Лоу (0) +5, / (0) +... + у" (Өш), 
Ка | T(n4 1) T (7) "PG, 
p @ T(n+r) Jo 


r being any positive quantity. 
[If + > 1 both integrals generally = zo .] 
[WorsTENHOLME, Educ. Times. } 


30. Prove by changing the order of integration or otherwise 


that z і 
а f ER [H5 e Uto) - 10). 


о Vy-£ 
(Maru. Trrpos, 1875.] 
31. Show that 


g^ nz" (m+ 1)2" 
de — ж n+l (m+1)(2n+1) (2n+1)(3n+1) 
huie ~ Wee ee ele ТС. 
(2n)2a" (2n + 1)?a" 
(3n + 1) (4n +1) (4n 4- 1) (5n 4+ 1) 
UCCA дает ем, 


[Lacnorx, Calc. Dif., vol. ii., p. 292.] 
Deduce expressions for log 1 + and tan-!z as continued fractions. 


32. Prove that 


2ur 


© 3 ро 

HC +5) - rir) T'(zo) Г (27), where ө=е%. 

[Sr. Jouw's, 1891.] 

33. Evaluate the modulus of Г(1+\/- Та).  [Swrrm's Prize, 1875.] 


34. Show that for very large integral values of m, Г(п+1) is 
very nearly the geometric mean between Г (п) and T (n + 1). 
[Охғовр, 1892.] 


35. If b be a large whole number, show that, provided z > — 1, 
T'(b 4-1) 
Г(=+1)' 
[ре Morean, Diff, Calc., p. 585.] 


(z - 1)(z 2)... (5+6) = т very nearly. 
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36. Writing ф(х) = е. z!/A/2x z**!, prove by the aid of Wallis’ 
theorem that ф (22) = [ф(2)]? when z is large. 
Then show that for any value of z, 


(a) p(z) recette) 


$(r-1) 
p(z) MAE E (n-1) (-1)" 
(0 lg 1) Tat 123 408 tone) z^ 577 
ф (ж) s $(2) 1 
(с) grin а (9) юв 1) Iaer 
$57 et ee TEAN» la 


(0) $22) -f4 С++ t CT 


where 6,, Ө,, 0,, ... are numbers between 0 and т. 
8 


A AATF @<#<+т» 


and finally deduce Stirling’s theorem, 
1.2.3... c= 2m e-*3* (1 + ex), 


where є; denotes a positive quantity which vanishes when 2 = о. 
[SERRET, Calc. Intég., p. 207.] 


37. Show that, if x be a whole number, 
logT (z+ 1) = 1 log 2r -x + (z +4) loga 


+E | (etm +5) log(1 *zis)- T 


[GuDERMANN.] 


38. Show that 
1.2.3...2> J2mz a*e-* and mate 1 


when 2 is large. [SERRET, Calc. Intég., p. 213.] 
39. Writing 
Dy mmrl ! 
$m Lt, EEE, qu, and ы-у УГ, 
(mn)! m ? 
prove that 


Uy, =U", „=“ Jas, $(n) 2 nb (27) F. 


Hence deduce Gauss' theorem, 
n 


weT() T (2+1). ‚г(«+®—\- (21) * n3T (nz). 


[ЗЕВВЕТ, Cale. Intégral, p. 190.] 
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40. Prove that 
e 1 ix i ve }® T 
> lie | T(n) —, (log vy? dv. 
Jum д, Morean, Diff. C., p. 594.] 


41. Prove that 


d * dt 
i, log TG) - logo | 16-01 «rp 


1 €x T] Wee 
pire crc ДЕШ) 
where C' is a certain constant. [Матн. Твтроз, Pt. II., 1915.] 


42. If the binomial expansion for a positive index be written 
(а de by nA 22 () qus. 
show that x( Bo -r4l, r+1)=1. 


Prove also that 

9r _, (1? (2)? (3)? (4)? 
| ТТГ eo 

43. Show that (1000)! lies between 
4:02387 x 107567 and  4:02388 x 102567, 
and is a number with 2568 figures in the ordinary system of 
numeration, its logarithm being 2567:6046442.... 
(Cournot, Théorie des Fonctions, vol. ii., p. 472.] 


44. Show that if 


A. В, 
log Г (£+ 1) =ю5\/2т + (z 4- 1) log z - 2+1—55 сщ р 348+. 
E m п—1 2n—1 e еее НЕШ 
= 3) Gn — D) nai + ( ч (22+ 2)” 
then _ R= | е7 a2" f2"+2 (ба) da, 
0 


а 
where f(a)= л 
[LrovviLLe, Journal de Mathématiques, Tom. iv., р. 317.] 
If Аһ be the maximum numerical value of f?"*?(a) between the 
limits «—0, а=, show that 
Е Ань 1 
(2n -- 2)! < (9n +1)(2n + 2) =" 
and examine the nature of the approximation attained by the 


omission of all the terms which contain Bernoulli's coefficients. 
[LrovviLLE, J. de M.; also Cournot, Théorie des Fonctions, p. 474.] 


and @ is a positive proper fraction. 
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45. Starting with 
" _„ e^ —-e-*P]dB 
вг) |7 СЕ вс “ГР В 


-Í (Р +Qe-) dB, say, 
0 
and putting A for the two terms with negative indices in the 


development of Q in ascending powers of В, namely pt 3p let 


в) [| (P+ Ве-=в) 4В and == (Q – К)е- ғар. 
0 D 

Then show that и 

(1) &() - 1 log 5. (2) Р) = 3 1087. 


(3) Р(х) - F(2)-1-z-(z-1i)logz. (4) P(x) =е-та7- 2m ен), 
(5) That when z is large е=(2) differs but little from unity. 
(6) log I'(z4- 1) 2 1log 2z + (2+1) n 2 


1X ceca у=, 
«p (= cy oo a5)? rir and 
(7) Deduce the equation, 


B 
log I'(z-- 1) = $log(27) + (2+ $)loggr-2* 155; 73-4 3+ 


Bao +(- 1 —— E Damn 1 us EN 
(2n — 1) 2n жа (2n + 1)(2n+2) дат, 


DAUO-IL [BeRTRAND, Calc. Intégral, p. 265.] 


*t(-1) 


46. Show that 


а) vf (pee 


(2) legT'(z 4 1)- | ^g as] ap. 


[Торномтев, Int. Calc., p. 392.] 
47. If A, be the acute angle whose tangent is the n power of 


the reciprocal of the 7% of the prime numbers 2, 3, 5, ..., show that 
24 2A 24, cos 2A =2 Bon 1 (2n)! За 
COs 2 1 COS & 2 COS = gC k^ 47-10 2 В,2 (4n)! , 
where B, is the n number of Bernoulli. [Млтн. Твгроз, 1897.] 
1 
48. If I= = show that 


T(1)-452333]3, (3) = rttr, 
T(3)225253-37-5 r(5)- тё 3-4-1, 
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1 ; 1 
49. If I= === and J= ; Ul ‚ show that 
(= 218244515 "ig iT iji TG)- ntis "ig ИЛЬ 
Г(у3)= тї” /:5ig ig ig t 1% J, TGs) = at gig is t rt Jt, 
SA s 2т МЕ Sor . 4m 
where $,-sin To’ S, = sin <= 10? 5, = sin ту, S, = sin т, 


and write down the values of D(45), P(j;) Г), Г(%), in 
similar form. 
Trt 
360° 
[Oxron» I. P., 1914.] 


50. Show that | z|log(1 +e) - e] dx = 
0 


51. Prove that the volume in the positive octant bounded by the 
planes 2 = 0, у = 0, z=h and the surface 2/с = a?" [am + ут" is equal to 


һу" irG)r 
“нага 


[Матн. TniP., Part П., 1913.] 
52. Prove that e ә {$(x)} -zf e~} p(x + 2y /h) du, 
and apply the result to prove that if 1 + 44k be positive, 


qi kg? 


e di { ye- Мыл} dl EE x Lg СТАЕ. 
(1 + 4hk)? 
= [Матн. Твте., 1870 (WOLSTENHOLME). ] 


53. When п is a positive integer, we have evidently 
1.3.3.5; 30425. 22.,20.1.2 (8 - Dj 


prove that this equation, when expressed by means of the function 
Г, is true for any positive value of n. [Sm G. б. Sroxzs, S. P., 1870.] 


94. Prove that the limiting value of 


i L (2n + 1)^ 
Colitis d v РР 
when n is indefinitely inereased, is log 2. [R. P.] 
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But there is this difference between the functions 


AC ANS UIT: " 
= G4-1)(z4-2)...(z4- и)“ 
that though they coincide in value for all positive values of z, 
the former becomes infinite at the values z— —1, z— —2, 
x= — 3, etc., but has finite values for other negative values of 


z, whilst the definite integral is permanently infinite for all 
negative values of 2-1. 


Live and 7 
о 


888. That the factor form has finite values, when p becomes infinitely 
large, for negative values of х between the asymptotes may be made 
clear by taking а case. Take r= — $. 


Lu2;2... 


Then Lt,-« 


2.4.6... р 1 
1:1.3.5...( 2n — 3) n. 


а ln 92, 42.62... (2)? (901) 
1.2.3.4... (2p 302и 2) (8 1)(8и) i 


(J Qua pe ^Y (24-1) 


= – 49 ———® : 
м Au (25) e7 sje m 
_ 27E Qrp 2ң— z s. - Ns. 
9 mp и? 
5 23 
Similarly at х= -z the corresponding limit is I3 Мт, 

at z= —7 the corresponding limit ів e UM - 
PAS Оа R34 


and so on. 
These mid-ordinates, half way between the successive asymptotes, thus 
form a regular descending series 


S4 ug ЦА. $9. Бє 2 
-ivm 1:59 РМ Pes etc. 
889. It is worth noticing that II(z, и) may be written as 


| uy aT: 
П (2, DI CESITCESJTCES NN CEU D 


OG GENCE» y 
(1+2) (1-3) (14+§)---(1+2) uFi 
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